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FOREWORD 


While Titan wasunder Contract to the NASA Dryden Flight 
Research Center (Contract NAS2-11990) , Titan was also under 
contract to the NASA Ames Research Center, Moffett Field 
(Contract NAS2-11824) . Both contracts had as ultimate goals to 
improve the Adaptive Maneuvering Logic Air-to-Air Combat 
Computer program. 

The emphasis in the Moffett Field program was to improve 
the guidance laws, regardless of required execution time on a 
computer. In contrast, the Dryden effort was to provide a robust 
decision logic, guaranteed to work in real time. The logic 
developed for Dryden should eventually drive an actual aircraft 
in real flight. 

During the course of this work, it would have been 
unproductive to keep book which of the AML improvements should 
be credited to the Moffett Field contract and which ones to the 
Dryden contract. This final report on contract NAS2-11990 is 
therefore essentially the same as the final report on contract 
NAS2-11824 (Simulation of Modern Air-to-Air Combat) . The present 
report has some material added in section 3 and a substantially 
enlarged section 5. It also contains an Appendix with the 
Fortran listing of the subroutines implementing the ’’Basic 
Fighter Maneuvers". 
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SIMULATION OF MODERN AIR COMBAT 


By George H. Burgin and Laurent B. Sidor 
Titan Systems, Inc., La Jolla Ca. 

SUMMARY 


This final report on Contract NAS2-11824 is organized 
in seven sections plus a list of 25 references. 

Section 1 provides an overview of current topics in the 
simulation of air-to-air combat, touching on such subjects as 
weapons simulation, aircraft modeling and performance measurement. 

In section 2, the history of a set of computer programs, 
developed over the last 15 years is traced. These programs are 
generally known as "Adaptive Maneuvering Logic" (AML) programs. 
They exist in many versions: Air-to-air combat and missile 

evasion, real-time and non-real time versions. 

The air-to-air combat simulation exists in two basically 
different versions: The older version, the "trial-maneuver" 

version, is described in other NASA reports. The newer version, 
the "IF ®> THEN " version, is the subject of section 3 of this 
report . 


Section 4 summarizes some important aspects of 
dynamics modeling. The interrelationship between the 


aircraft 

tactical 


Vll 



decision process and the aircraft model is shown. For example, 
tactical performance can be significantly improved if the 
modelled aircraft can be controlled by a control system capable 
of orienting the aircraft's longitudinal axis into a desired 
direction (pointing control system) . 

Section 5 compares the performance of the trial maneuver 
logic with the IF => THEN logic and demonstrates how each logic 
may be improved by "playing" it against the other logic. 

To make the performance of the IF *> THEN logic less 
predictable, some basic fighter maneuvers were added to AML which 
are invoked, when appropriate, under the control of pseudo- 
random numbers. These maneuvers are described in section 6. 

Finally, section 7 provides some suggestions for 
continued work in developing advanced guidance law for air-to- 
air combat. 
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1. OVERVIEW OF AIR COMBAT SIMULATIONS AND METHODOLOGIES 


Simulation of air-to-air combat has become an 
indispensable tool for pilot training, for tactics development, 
for weapons systems evaluation and for a host of other 
applications. Air combat becomes more and more complex due to 
advances in electronic warfare. Air-to-air combat today begins a 
long time before the opponents have visual contact. Radar and 
other sensors provide critical input to the pilot at a range far 
beyond the visual range. It has therefore become common practice 
in the analysis and in the simulation of air-to-air combat to 
differentiate between a “Beyond Visual Range ( BVR) '• phase and a 
"Close In Combat (CIC)" phase. 

The present report is concerned primarily with simulating 
the CIC environment. Specifically, we will discuss in detail a 
series of computer programs generally known as "Adaptive 
Maneuvering Logic Program (AML) . These models and simulations 
were developed under NASA sponsorship with the initial goal to 
have an intelligently interactive, real time opponenet on NASA’s 
differential maneuvering simulator (DMS) at the Langley Research 
center and with the long-range goal to provide assistance to a 
pilot during air combat engagements. 

A measure of the complexity of modern air war may be 
obtained by reading the account of Israeli air operations over 
Lebanon in 1982 (Reference 2) . These operations involved air 
superiority fighters in strike escort and combat air patrol 
roles operating in concert with many other elements such as 
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SAMs, AWACS, ground-based radars and communication centers, 
stand-off jammers, and RPVs. Similarly complex operations are 
involved in the air defense of U.S. carrier task forces (see 
for example Reference 3) . 

Due to the complexity of such air operations, individual 
air simulations must focus on a particular, limited area of air 
combat. We will briefly review the current state-of-the-art in 
air combat simulations in order to put a perspective on the area 
considered by AML. Some of the key issues addressed by this 
report will be: 

- Number of aircraft involved in the simulation 

- Types and properties of weapons employed by the 
combatants 

- Degree of complexity of aircraft and weapons models 

- How random effects are simulated 

- Off-line simulation versus real-time simulation 

1.1 DETERMINISTIC AND STOCHASTIC MODELS 

The common point of departure for air combat simulations 
are various scenarios of Air Force and Navy Missions. In the 
final analysis, their common evaluation point relies on pilot 
opinion. In the design phase, a basic trade-off must be made 
between the accuracy in modelling individual elements and the 
size and execution time of the code. Figure 1.1 attempts to 
portray this trade-off. Engineering simulations which model in 
detail physical mechanisms (such as warhead fuzing) are limited 
to one or two units. At the other end of the spectrum are 
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LEVEL OF DETAIL IN SIMULATIONS 
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Figure 1-1. Level of Detail in Simulations 
(Adapted from Reference 1) 



campaign or force-on-force models with hundreds of simulated 
units. In these models, the representation of physical 
mechanisms in the simulation is done in terms of aggregated 
performance measures. The simulation of even a minimally 
representative number of opponents in the case of NATO vs Warsaw 
Pact scenarios (2 vs 4) leads to an explosion in the 
computational requirements . 

The performance of many aspects of weapons systems is 
expressed in terms of probabilistic quantities, for example 
radar probability of detection or kill probability of a missile 
warhead against a target type. The combination of these 
probabilities can be performed in one of two ways: (1) Expected 
value method and (2) Monte-Carlo method. In the expected value 
method, the probabilities are combined using the law of 
probabilities for the particular probability law obeyed by the 
simulated process. For example, if there are N independent 
interceptors, each with a probability PD of detecting a single 
bomber over a period of time, then it may be shown that the 
expected fraction of bombers FDB detected at the end of the 
period of time will be : 

FDB - 1. - exp( N *PD/ M) , where M is the total number of 
bombers. (Reference 15) 

In contrast, in the Monte-Carlo method, the outcome of a 
probabilistic event is assessed based on the draw of a random 
n umb er - For this reason, these are described as "discrete 
events". Repetitive trials must be performed to obtain averages, 
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a process which multiplies computational requirements. In 
addition, in Monte-Carlo simulations, the sheer volume of 
information makes it difficult to trace causative factors. For 
these practical reasons, Monte-Carlo models are popular up to 
the mid-range of Figure 1-1? for campaign models, only expected- 
value models are practical. 

1.2 OFFLINE AND REAL-TIME SIMULATIONS 

The simulations discussed above consist of "off-line" or 
"non-real time" simulations. Even these non-realtime 
simulations may require execution times which limit their 
economical use for studies and analyses. 

In a real-time simulation, two tasks have to be performed. 
First, the equations of motion for each participant must be 
solved satisfying the condition that the CPU-time to perform the 
calculations for one integration step will not exceed the 
allocated frame time. The second task requires the simulation 
of the decision process for each platform. The required CPU time 
to perform this second task also must fit into the allocated 
frame time. Typical frame times for real-time close-in air-to- 
air combat simulations with a human pilot in the loop are 
between 10 and 50 milliseconds. 

Table 1-1 illustrates the parameters involved at both ends 
of the spectrum in complexity in air combat simulations. The 
AML programs feature a high complexity simulation environment, 
moderate complexity in aircraft performance and tactics 
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representation, and low complexity in weapons and avionics. 

1.3 WEAPONS MODELS 

The armament considered in air combat simulations consists 
of guns, guided missiles and, recently, lasers. Because of the 
research nature of these simulations, a significant amount of 
effort has been spent on simulating air-to-air lasers, while 
this weapon has yet to see operational use. They will not be 
discussed further here. 

The basic requirement to achieve a gun firing position is 
to point the nose of the aircraft at the target. Steerable guns 
would alter this requirement, but presently there are none 
operational on fighters anywhere in the world. In general, to 
achieve a kill will require several hits on the target. For this 
reason, an off-tail position is preferable (Figure 1-2) . In AML, 
the conditions to achieve a gun firing position are a line-of- 
sight angle less than 10 degrees and an angle-off tail less 
than 60 degrees and a range less than 3000 feet. (These 
quantities are defined in the paragraph "Performance measures" 
below. ) Some models provide the option of integrating the 
trajectory of an individual bullet. The point of impact is 
calculated so that the effect of the hit can be accurately 
estimated using a vulnerable area approach. 

The requirements that must be satisfied for a missile 
launch are customarily summarized in terms of a "firing 
envelope" or "launch-acceptable region (LAR)". A representative 
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Cone of Vulnerability 



Above: Most successful gun 
attacks are maije from astern at 30* 
angle-off or less, the lethal cone. 
Vulnerability cone is 45* angle-off 
and 1,640yd (1,500m | range. 


Figure 1-2: Gun Vulnerability Cone 

(from Reference 24) 
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envelope, with the target at the center, is shown in Figure 1-3 
for a typical radar doppler-homing missile (Reference 10) . These 
envelopes are used by most m on n air combat models. They have 
also been used in various studies using AML (Reference 13) . 

For non-maneuvering targets, an envelope has a maximum 
range with a roughly elliptical shape which reflects the aero- 
propulsive limit of a typical missile. The maximum range varies 
strongly as a function of altitude and target speed. Figure 1-3 
b also illustrates the seeker limit, which in the illustrated 
case is smaller than the maximum range of the missile. The 
seeker limit is dependent on the target's radar reflectivity 
characteristics (a function of the target aspect as seen from 
the firing position.) Figure 1-3 also indicates an inner zone 
(minimum range or "dead zone"). It will be noted that the head- 
on maximum range is much greater than the off-tail range — 
typically four to five times. The greater area means that there 
are more engagement opportunities in the forward target quarter. 
But it should also be remembered that the target's sensors are 
effective only in its forward quarter. 

The "maneuvering envelope", as illustrated in figure l-3b, 
represents the effect on the intercept capability of the missile 
when the target begins a level left-hand turn just as the 
missile is launched. The envelope typically assumes a shape that 
is distorted in the direction of the turn. The maximum range 
expands in the direction of the turn, as the target is flying 
towards the missile. It contracts in the direction opposite the 
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MAXIMUM 

RANGE 


a) Non-maneuvering Target 



b) Maneuvering -Target 


Figure 1-3: Typical Missile Envelopes for Maneuvering 

and Non-maneuvering Targets (from Reference 10) 
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turn as the target is flying away from the missile. These 
effects are proportional to the number of G's pulled by the 
target. In spite of the magnitude of the effect due to 
maneuvering on the shape of the envelope, none of the models 
listed in the references appears to modify the decision to fire 
based on this effect. 

After the missile has been fired, the damage to the target 
must be assessed. Detailed simulations simulate the fly-out 
trajectory to the target, compute the miss distance and 
resulting survivability of the target. Less detailed simulations 
simplify this problem by computing a time-of-flight and 
survibability of the target based on a probability of kill and 
Monte Carlo draw. 


1 . 4 AIRCRAFT MODELS 

The simplest aircraft model used in air combat simulations 
consists of a point mass to which are applied the lift, thrust 
and drag forces. This provides a starting point, for instance to 
compute the endurance of an aircraft in the simulated 
engagement. This type of aircraft model is limited to 
"instantaneous turns", and cannot represent the attitude and 
turn capability of a fighter. Yet this limitation is often not 
recognized until realistic graphics are available, or the 
simulation is run in a flight simulator. 

A full 6 degrees of freedom (6 DOF) model is required to 
simulate realistically the behavior of a fighter aircraft under 
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the high-Gs and very large angles of attack encountered in air 
combat. All the lift, drag and thrust characteristics as well as 
moments should be represented. In particular, roll performance 
is of primary importance in fighter aircraft tactics and would 
alone justify the use of a 6-DOF model (see for example 
Reference 18 . ) 

The model currently used in AML is described in Reference 
14. It is a "performance model", in which 6-DOF dynamics have 
been preserved, but in which the calculation of aerodynamic 
moments and control and stability derivatives has been omitted 
to meet execution time requirements on minicomputer-based flight 
simulators. 

1.5 DECISION-MAKING PROCESS 

The objective of the decision-making process is to derive 
maneuvers which will bring one's own weapons to bear on the 
target while at the same time minimizing exposure to the other 
side's weapons. It is essentially a representation of the action 
of the pilot during combat. In simulations involving multiple 
aircraft, the decision-making process also involves pairing 
groups of opponents. 

The real-life approach to the solution of the problem of 
steering an individual aircraft relative to an opposing, 
dissimilar aircraft is known as the "Basic Fighter Maneuvers" 
(BFMs) . Examples of training manuals describing BFMs for 
particular aircraft may be found in References 11 and 12 . 
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Reference 10 is a more general treatment of this field. The 
objective of BFMs is twofold: (l) Gain and maintain a positional 
advantage with respect to the enemy allowing employment of 
armament, and (2) Gain and maintain sufficient energy to have 
maneuvering potential. BFMs are not exact maneuvers, but rather 
combinations of the three elementary actions that an aircraft is 
capable of — roll, turn, and accelerate/decelerate — used to 
gain advantage in a particular situation and against a 
particular opposing aircraft type. Well-known examples of BFMs 
are: the Immelman, the lead/lag turn, the Lufbery, the high- 
speed yo-yo. 

In spite of the admittedly inexact nature of BFMs, they 
nevertheless constitute a sourcebook of possible maneuvers which 
has been used as the basis for the decision logic of models such 
as PACAM (Reference 8) , AASPEM (Reference 7) and TACBRAWLER 
(Reference 9) . As an example of this approach, a partial list of 
such maneuvers available in the AASPEM model includes: 
chandelle, split-S, high-speed yo-yo, barrel roll. The 
decision logic for selecting a maneuver is based' for the most 
part on user-specified geometry rules. There is an amount of 
guesswork involved in specifying these maneuvers. For example, 
the user must insure that the energy state of each aircraft is 
sufficient to complete the specified maneuver. otherwise, 
unrealistic and unacceptable maneuvers may result. 

The specification of these maneuvers depends on the current 
phase of the engagement. For example, AASPEM considers seven 
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phases: 


- Neutral: no threat detected 

- Late set-up: setting-up phase near completion 

- Early set-up: setting-up phase near completion 

- Pre-attack: final set-up and preparation for attack 

- Attack: attacking threat 

- Post-attack: initial attack complete 

- Disengage: engagement complete 

For each of these phases, AASPEM requires specifying 
positional tactics, information-gathering tactics and 
information-denial tactics. 

This approach suffers from the disadvantage which was noted 
in the original AML report (Reference 5) , and is echoed in some 
training manuals (Reference 11) that fighter pilots learn these 
basic fighter maneuvers in training, but they rarely complete 
them in a dogfight because of the continuous interaction and 
changes in the relative situation. 

Another type of approach consists of programs which apply 
such disciplines as optimal control theory, and the theory of 
differential games to obtain control laws. Such approaches work 
best for idealized situations (e.g. co-altitude, analytic lift 
curves, etc. . . ) . 

The trial-maneuver approach was introduced by AML to remedy 
the problems with these approaches. The AML technique determines 
the next tactical maneuver as it contributes to the goals of the 
pilot. It uses the concept of a situation matrix describing the 
tactical decision options in terms of various values assigned to 
each cell. The maneuver selected is the one which maximizes this 
value (References 5 and 6) . 


1-14 


The trial-maneuver approach, originally published by Burgin 
et al (References 5 and 6) proved to be quite successful for a 
real-time decision logic. It does require, however, considerable 
computer resources. Pedotti and Hignard (Reference 22) 
plagiarized the above mentioned work. They used almost an 
identical set of trial maneuvers and had a real time version of 
their " Logique Adaptive de Manoeuvre Aerienne" running on an 
UNIVAC 1100/82 mainframe computer. 

Austin et al ( Reference 23) used a very similar trial 
maneuver technique in the simulation of air-to-air combat 
between two helicopters. This program is operational in real 
time on the NASA AMES Vertical, Motion Simulator. 

The trial-maneuver approach — as the name implies — 
involves searching over a series of flight paths. The 
computational requirements were found to exceed the capacity of 
VAX 11/780-class mini-computers for real-time applications. To 
remedy this situation, a different approach was devised: the 
rule-based AML (RB/AML) . The rule-based AML uses a combination 
of production rules (i.e. IF ... THEN statements) and guidance 
laws as an alternative to the trial maneuvers. These rules will 
be discussed in greater detail in Sections 3 . 

In m-on-n simulations, the decision process must in 
addition pair various groups of opposing aircraft. The doctrines 
found in the tactics manuals are the welded -wing, free-engaged, 
and the double attack system. These tactics have been emulated 
in air combat models such as PACAM and AASPEM. In the welded- 
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wing doctrine, the wingman attempts to maintain a loose 
formation with his leader. He does not make independent maneuver 
decisions, but nevertheless he fires his weapons on his own 
initiative when such opportunities arise. In the doctrine of 
free-engaged tactics, the two fighters exchange the roles of 
leader and wingman as the tactical situation requires. 

1.6 PERFORMANCE MEASURES 

Three levels of performance measures can be found: 

(1) Individual aircraft performance, e.g. turn rate or 
energy maneuverability as a function of Mach number and 
altitude. The ability to change state is a recently introduced 
performance measure in this category (Reference 18.) 

(2) Differential aircraft performance measure, e.g. 

the difference in turn rate. These are commonly used in training 
manuals. The implicit assumption is that both opponents enter 
the combat arena under the same initial conditions. 

(3) Tactical performance measures, which are made possible 
only through air combat simulations of the type analyzed in this 
report. 

The relative position of two opposing aircraft, "A" and "B" 
is conventionally described in terms of the deviation angle 
lambda and angle-off epsilon. These have been illustrated, from 
the point of view of "B", in Figure 1-4 , where they are 
indicated as lambda(B) and epsilon (B) . The deviation angle 
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lambda (B) is the angle between "B M, s velocity vector and the 
line of sight from »B» to ''A". For this reason, it is sometimes 
referred to as the "line-of-sight angle". This deviation angle 
is an indication for "B" of where "A" is: lambdab - 0 degrees 
means "A" is directly in front of "B"; lambdab - 180 degrees 
means "A" is directly behind "B". 

The angle-off epsilon (B) is measured between the line-of- 
sight vector from "A" to "B" and "A"'s velocity vector. It tells 
"B" where "A" is going relative to "B": epsilon (B) - 180 degrees 
means "A" is coming directly at "B"; epsilon (B) - 0 degrees is 
going away from "B". Alternate names for angle-off are: angle- 
off-tail and aspect angle (Reference 11, page 2-2) . 

Similar angles can be defined for "A". Inspection of figure 
1-4 shows that lambda (B) * 180 deg - epsilon (A) and lambda (A) = 
180 deg - epsilon(B). 

The line-of-sight angle and the angle-off are 
fundamentally important in air-combat; both for the tactical 
decision process as well as for the assessment of the current 
situation. A few clarifying remarks are therefore in order. 

First note that the AML program carefully differentiates 
between line-of-sight angle and deviation angle. In the follow- 
ing discussion, we reference all the angles to aircraft "B", in 
other words, when we say, line-of-sight angle, we mean aircraft 
"B" 's line-of-sight angle. By AML's definition, the line-of- 
sight angle is the angle between the vector from "B"'s eg to 
"A"'s eg (the line-of-sight vector) and "B'"s body x-axis . The 
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"B"'s Perspective: 

o The deviation angle XB tells "B" where "A" is. 

o The angle-off £B tells "B" in which direction "A" is 
going with respect to "B". 


Figure 1-4. Definition of Deviation Angle and 
Angle-off ("B"'s View) 
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deviation angle, on the other hand, is defined as the angle 
between the LOS vector and "b ,M s velocity vector . Line-of-sight 
angle and deviation angle therefore are only identical if there 
exists no sideslip and no angle of attack. 

From a tactical point of view, both the deviation angle and 
the angle-off are important. For gun-firing, the line-of-sight 
is of primary importance, because the guns are mounted such that 
they point in the direction of the aircraft's longitudinal axis. 
For missile firing, both the line-of sight angle and the 
deviation angle are important, the missile is mounted parallel 
to the aircraft's longitudinal axis, the initial missile 
velocity, however, is determined by the aircraft's velocity 
vector. 

One last point: The line-of-sight vector can be changed by 
the pilot much more rapidly than the deviation angle, because 
modern fighter airplanes allow very rapid changes of angle of 
attack of the order of 10 to 20 degrees. This translates 
directly into a line-of-sight angle change of the same 
magnitude. The velocity vector however can not be changed that 
rapidly. 

Although these definitions (or equivalent definitions) are 
widely in use air training manuals as well as air combat 
simulations, it should be noted that there are ambiguities 
arising from the fact that the values of the line-of-sight angle 
and the angle off are between 0 and 180 degrees and always 
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positive. 


For example, if one expresses the situation in the 
lambda (B) -epsilon (B) plane, both situation 1-5 a and 1-5 b will 
be represented by the same point in that plane, namely lambda (B) 
=*90 degrees, epsilon (B) =90 degrees. 

Similarly, the two situations 1-5 c and 1-5 d fall into the 
same point in the lambda-epsilon plane, lamda(B) =90 degrees, 
epsilon (B) = 45 degrees. 

Since 1-5 a represents a tactically different situation 
from 1-5 b, these ambiguities should be removed if one wants to 
base the tactical decision on the two angles lambda and epsilon. 
One possibility would be to introduce, in addition to these two 
angles, also the line of sight angle rate. Assuming equal 
velocity for the two aircraft, the line of sight angle rate 
would remain zero for situation 1-5 a, but a large rate would 
result in situation 1-5 b. Similar observations can be made 
between situations 1-5 c and 1-5 d. 

Notwithstanding these limitations, it has been found useful 
to introduce a performance index which combines these two angles 
into a single measure: 

PI (B) = 50* ( 1 - lambda (B)/180) + 50*(1 - epsilon(B)/180) 

PI (A) = 50* (1 - lambda (A)/ 180) + 50* (1 - epsilon(A)/180) 

Values of PIB and PIA are illustrated in table 1-2. 

For example, during its AML tests (Reference 4) , 
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1-5 c 1-5 d 


Figure 1-5 Ambiguities in the Line-of-Sight/Angle-Off 
Representation 
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NASA/Langley considered that a pilot (say "B") enters a gun zone 
if lambda (B) does not exceed 10 degrees and if epsilon (B) does 
not exceed 60 degrees (and the range is less than 3000 feet) . 
This condition corresponds to a performance index PI(B) of 80 
or better. (Correspondingly, PI (A) would be 20 or less.) It 
should also be noted that this value of the PI is a necessary 
but not sufficient condition for a gun-firing position. 

The integrated performance index is the time-averaged value 
of the instantaneous values of the performance index: 

IPI = ( PI1*DT + PI2*DT + + PIN*DT ) / T 

Offensive time 

The offensive time is defined as the accumulated time 
during which the opponent was in front of the wing line of the 
reference aircraft. This was one of the figures-of -merit used 
during the original AML test runs at Langley (Reference 4) . 

The offensive time with advantage is the accumulated time 
during which the opponent was in front of the reference 
aircraft's wing line and the reference aircraft was behind the 
opponent's wing line. AML also used a more restrictive 
definition of the offensive time, consisting of the accumulated 
time during which the reference aircraft's deviation angle was 
less than 60 degrees and its angle off less than 60 degrees. 

Time to first kill 

Other performance measures account for the weapon 
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probability of kill. The time to achieve the first kill is an 
attractive measure of this kind. However, its drawback is that 
it does not properly reflect the future impact of the aircraft 
state at the time the first kill is achieved and its subsequent 
capability to engage more targets. An initial firing position 
may be achieved by turning at maximum instantaneous load factor 
in order to gain an angular advantage. However, this will result 
in the aircraft losing rapidly energy and thus position itself 
unfavorably for a subsequent engagement. In typical air-to-air 
scenarios, it is precisely the purpose of the leader/wingman 
team concept to take advantage of such situations. 

Accumulated probability of kill 

A commonly used measure of military effectiveness is the 
loss exchange ratio, defined as the ratio of the number of enemy 
killed divided by own losses. In a one-on-one duel in which 
multiple weapons are exchanged, this reduces to the ratio of the 
accumulated probabilities of kill. Neuman and Erzberger used 
this measure (Reference 13) as an alternative to the measures 
of effectiveness previously discussed. 

The common procedure to calculate the exchange ratio is to 
use the Monte-Carlo method. An alternate method was used in 
Reference 13 in which the engagement continues independently of 
the outcome and these trajectories are recorded. A post- 
processing program uses these trajectories to identify firing 
opportunities and to compute the accumulated probability of 
kill. This method was used to avoid the problem often 
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encountered in air combat simulations that small changes in 
initial conditions or in the flight path somewhere in the 
engagement propagate into large differences in outcome. However, 
this method is limited to lvl, since in m on n there are 
cooperative effects which depend on the sizes of the forces. 
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2 . DEVELOPMENT HISTORY OF THE AML PROGRAMS 

2.1 NASA LANGLEY DMS PROGRAM (AML 75) 

Development of the AML program started in 1969 under the 
sponsorhip of the NASA Langley Research Center. The original AML 
program was developed to operate on NASA Langley Research 
Center's Differential Maneuvering Simulator (DMS). It is shown 
in figure 2-1 as the root of the AML program family. It was 
designed to be an interactive air-combat opponent operating in 

real time. This original version of the AML program is 
documented extensively (for example, References 5 and 6) . 

2.2 THE DMS CONTROL MODEL (AML 76) 

In the original version operating on the DMS, AML would 
"drive" the displayed aircraft by providing body rotational 
commands p, q, and r to the DMS display program. AML calculated 
the values for p, q, and r such that the displayed aircraft 
would achieve an attitude compatible with the following 
conditions: 

- resultant force vector ( aerodynamic forces, propulsive 

forces and gravity force) must 
lie in the desired maneuver 
plane. 
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Air-tQ-Air .Combat 


Missile Evasion 



Figure 2-1 . Development History of the AML Air-to-Air and Missile-Evasion Programs 
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- the angle of attack is such that the desired lift is 
produced 

- sideslip is zero 

The AML program then filtered these commanded body rotational 
rates in order to achieve a smooth motion of the displayed 
aircraft. However, when a pilot lost an engagement with the AML, 
he had a tendency to claim that the AML driven aircraft would 
perform flight maneuvers which were outside the performance 
envelope of the real aircraft. To counter this argument, a 
control system was developed which would actually move the 
simulated aircraft's control surfaces in exactly the same manner 
as the pilot did it with the stick and the rudder pedals. These 
commanded control surface deflections were then fed into the 
identical set of equations of motion as the were used to drive 
the human piloted aircraft. The development of this control 
system is described in reference 20. A thorough comparison 
between the performance of the original AML (called the 

performace-model AML) and the AML with a control-system is 
contained in reference 4. 

2.3 A ONE-VERSUS-TWO VERSION OF AML (AML 1V2) . 

The Human Research Laboratories of the Air Force (AFHRL) 
sponsored subsequently an extension of the one-versus-one AML 
version to a one-versus-two version Here AML represents the 
single aircraft opposing two bogeys. This lead to a batch- 
version of AML which handles the one-versus-two situation based 
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on a set of value functions taking into account the relative 
situation between all three combatants. This version also 
replaced some of the binary value functions of the original AML 
by continuous functions, thus giving' a better resolution between 
trial maneuvers and avoiding ambiguity in the scoring of 
different maneuvers. 

2.4 AML WITH REVISED EQUATIONS OF MOTION (AML 84) 

The original AML (AML 75) had a number of known 
deficiencies. The most serious one was an abnormal behavior of 
the AML aircraft when it approached 90 degrees in a vertical or 
in a near vertical turn. This anomaly was not due to the 
singularity of the Euler angles at theta = 90 degrees ( AML uses 
quaternions for the attitude integration and consequently there 
is no singularity at any attitude) . The problem rather had to do 
with the decision logic and it may be explained in somewhat 
simplified form as follows: Most maneuvers in AML are executed 
in "maneuver-planes". A maneuver-plane always passes through the 
aircraft's velocity vector. Certain other parts of the decision 
logic are based on the line-of-sight angle, which is, in part, 
determined by the direction of the aircraft's body-x axis. In a 
vertical loop, under high angle-of-attack conditions, it will 
happen that the body axis has already exceeded the 90 degrees 
pitch angle, but the velocity vector's pitch angle is still 
below 90 degrees and still increasing. Under such a situation, 
it can occur that the AML reverses its maneuver command 
inappropriately. Specifically, it will command a maneuver plane 
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rotation angle of 180 degrees (or close to 180 degrees) when in 
reality, the maneuver-plane rotation angle for the intended 
maneuver should be zero (or close to zero) . This reversal can 
take place during several subsequent decisions. The result in a 
flight simulator is that it looks as if AML wouldn't know what 
to do. The long range effect is a hammer-head stall of the AML 
aircraft. Section 4 of this report explains briefly how this 
problem was solved. 

2.5 NORTHROP AEROSCIENCES LABORATORY AML (INTERACTIVE 
TARGET) . 

The Aerosciences Laboratory of the Aircraft Division of the 
Northrop Corporation, which had an early version of AML 
installed on ther moving base simulator, was interested in an 
AML implementation with the new equations of motion, which 
eliminated completely the "over the top" problem explained in 
the previous paragraph. However, the computer hosting AML was a 
Harris Slash 4 minicomputer whose computational capability was 
inadequate to support AML in real time, not even with a frame 
time as large as 50 milliseconds. To reduce execution time, we 
abandoned the concept of trial maneuvers and of selecting the 
most promising of these trial maneuvers. Instead, we developed a 
logic which resembled closely the production rules of the then 
popular expert systems. This not only allowed us to perform a 
tactical decision well within the allocated frame time, but it 
also gave AML the flavor of an AI program. 


2-5 



2.6 GENERAL ELECTRIC'S AML VERSION ON THE SIMULATOR FOR 
ADVANCED AIRMANSHIP (AML 86) 

At the time of this report, this is the most advanced real 
time version of AML. It uses basically the same decision logic 
as the Northrop version, with some added improvements for low 
speed, low energy avoidence. The host machine is an SEL 32/97 
computer and the visual display is a General Electric Compuscene 
III computer generated image. AML86 has a number of additional 
features, such as minimum allowable altitude for the AML 
aircraft, a choice between three different aircraft (F-4, F-5 or 
F-15) , a large number of selectable "canned" maneuvers for the 
AML aircraft and most interesting, a selectable skill-level for 
the AML aircraft. The skill level of the AML aircraft can be 
selected to be "ACE", "AVERAGE", or "GRAPE". 

2.7 AML 87 / EXPERT IVAN 

Presently, there is an in-house effort going on at Titan 
with the two objectives of: 

1. - Expand AML's decision logic to BVR 

2. - Expand AML's decision logic to handle multiple 

aircraft on both sides. 

AML 87 is strictly a production rule based system, the rule-base 
being built by Navy fighter pilots with current experience in 
air-to-air combat in F-14's and F-18's. 
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2.8 MISSILE EVASION AML (AML/SAM) 


The initial success of the AML program as an "iron pilot" 
in the DMS created confidence that the AML decision logic could 
be changed to "fly" AML such that it would avoid a surface-to- 
air missile. This work was initially sponsored by the Tactical 
Fighter Weapons Center at Nellis AFB, Nevada. The obvious 
required change to AML was to replace the value functions (which 
favored achieving a six-o'clock position with respect to the 
opponent) to functions which favored achieving a large distance 
between the missile and the AML aircraft. Obviously, the type of 
trial maneuvers also had to be changed. Less obvious is the fact 
that in case of missile evasion, a short term maneuver 
optimization (as it is performed in the air-to-air combat 
version) will not generate maneuvers with accepable miss- 
distance. It is necessary to carry out the optimization from the 
decision time all the way to the impact (or the point of closest 
approach of the missile) . The decision logic of the AML program 
was modified to implement these requirements and very successful 
evasive maneuvers against surface-to-air missiles, such as the 
SA6 were generated by AML. The program ran in non-real time on a 
CDC Cyber computer. 

2.9 PILOT'S ASSOCIATE D1 AML PROGRAM ( SAML Dl) 

The Aircraft Division of the Northrop Corporation 
participated in the demonstration phase of the Pilot's Associate 
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D1 program. The AML/SAM program was modified to work in real 
time on a flight simulator. AML determined suitable evasive 
maneuvers for the aircraft. These maneuvers were generated based 
on a set of production rules. The AML generated maneuvers were 
either used to provide cues displayed to the pilot on the 
heads' up display or they were fed directly into a flight control 
system. AML successfully avoided, at very low altitude, two SAMS 
simultaneously. For further details, see reference 21. 



3. THE BASIC AML IF => THEN LOGIC 


3 . 1 TERMINOLOGY 

The purpose of this section is to give an overview of the 
tactics currently implemented in the IF => THEN version of AML 
which is in use at various flight simulation facilities. To 
avoid confusion, we will first clarify some terminology. In a 
real-time, one-on-one environment on a flight simulator, the AML 
driven aircraft is called A/C "B", or for short, AMLB. When the 
program operates in a batch environment, the opponent of AMLB is 
AMLA. Each of AMLA and AMLB can implement either the "trial- 
maneuver-logic" or the IF *> THEN logic. The following 
discussion assumes that the , 'B H aircraft is driven by the IF => 
THEN logic. In the rest of the report, AMLA is a "trial-maneuver 
logic " AML. 

3.2 COMMONALITY WITH THE TRIAL-MANEUVER LOGIC 
3.2.1 Timing Considerations. 

In the IF => THEN logic as well as in the trial maneuver 
logic, two time-intervals are used for maneuver decisions. The 
first one, which is the smaller of the two is equal to the 
integration stepsize ( alternatively called frame-time or cycle- 
time) . The AML maneuver logic subroutine (TACTICB, see figure 3- 
3) is invoked every integration step. At each invocation, the 
AML logic unconditionally checks for the necessity of either 
initiating a dive recovery or to continue a dive recovery 
currently in progress. If no dive recovery requirements exist, 
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then the logic tests whether there is tine to perform a new 
tactical decision. This second tine-interval between tactical 
decisions is called decision-interval. For close-in, one-on-one 
air-to-air combat, it is typically between 0.5 and 1.5 seconds. 

(For missile evasion it is shorter and for a decision in a BVR 
situation it may be considerably longer) . 

3.2.2 Maneuver Plane Concept. 

A significant contribution to the success of the early 
versions of the AML program came from the concept of the 
maneuver-plane. Strictly speaking, one should not call this 
plane a maneuver plane, but rather a maneuver half-plane. (As it 
is properly called in reference 22). It is the half plane in 
which, ideally, the next segment of the AML aircraft velocity 
vector will lie. It extends through the AML driven aircraft's 
velocity vector towards the side of the cockpit. The maneuver 
plane provides (1) a convenient mechanism to specify AML 
maneuvers ( both in the trial maneuver and in the IF -> THEN 
version) and (2) a computationally efficient way for prediction 
of the aircraft's future position and attitude. In the IF »> 
Then logic, the maneuver-plane serves to specify the parameters 
for lead or lag pursuit maneuvers. The maneuver plane and its 
associated maneuver plane coordianate system are illustrated in 
figure 3-1. 

The crucial problem in both AML versions is to control the 
aircraft's body rotational rates in such a manner that: 
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Figure 3-1. Illustration of Maneuver-Plane Concept 



a) they are physically executable under the prevailing 
flight conditions and 

b) the aircraft's velocity vector remains in the 
specified maneuverplane. 


3.3 DECISION HIERARCHY 

3.3.1 Ground Avoidance 

The ground-avoidance logic is executed every integration step. 
This reflects the fact that ground-avoidance has higher priority 
than any other tactical decision. In both AML versions the 
decision on whether a ground-avoidance maneuver is required is 
based on a two dimensional table of the dive recovery angle. 
This angle is a function of airspeed and altitude. In the IF => 
THEN logic, a dive recovery maneuver leaves no choice, it is a 
roll to wings level followed by a maximum instantaneous g 
pullup. The throttle is controlled such that the aircraft is 
going to fly at corner-velocity. The dive-recovery maneuver may 
therefore succintly be described as a maximum g turn in a 
maneuver-plane whose rotation angle rho is zero. 

3.3.2 The Pointing Algorithm 

If dive recovery is not required, the program performs a 
test whether the aircraft should be controlled in such a way 
that its nose (i.e. its longitudinal axis) will point at the 
opponent or at a specified point in front of the opponent. This 
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is the only maneuver in the AML program (except of course 
the "canned" maneuvers in certain versions on flight simulators) 
where the maneuver is not based on a maneuver plane, but where 
directly body rotational rates which will bring the aircraft 
intothe desired attitude, are calculated. The pointing algorithm 
is described in more detail in section 4. 

3.3.3 The Lead/Lag Maneuver Logic . 

These maneuvers form the heart of the basic AML maneuver 
decision logic. They implement one of the basic rules of air 
combat: Point your nose towards the opponent. The refinment 
consists in the determination of the exact point in reference 
with the opponent towards which we want to point the aircraft 
(behind = lag, in front = lead or exactly at the opponent = pure 
pursuit) ; the other refinments being the rate of turn by whcih 
we want to achieve this goal ( in other words, the loadf actor) 
and finally how much thrust we will apply (throttle setting) . 
The decision on whether to fly lead, lag or pure pursuit is 
based on the values of the line of sight angle and the angle 
off, as illustrated in figure 3-2. 

Load Factor Selection The load factor is also selected 
as part of the LLG. This selection process, however, is 
primarily determined on the basis of airspeed considerations. A 
high load factor results in a high turn rate, which is desirable 
to achieve a firing position as quickly as possible. However, 
turns at the maximum load factor create a lot of drag which 
causes the airspeed to drop rapidly. This is actually desirable 
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Figure 3-2. Regions for Steering Laws 
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when the current airspeed is above corner velocity, the velocity 
which yields the highest turn rate. For this reason, when the 
current airspeed is above corner velocity, the maximum load 
factor is commanded. When the current airspeed is near or below 
corner velocity, the sustained load factor is commanded to avoid 
losing further energy. In B's forward sector (LOS < 60 degrees), 
an additional test is performed which compares the load factor 
described previously, which is airspeed-oriented, to the load 
factor corresponding to the desired flight path, i.e. the flight 
path which intercepts the reference point. This "intercept 
trajectory" load factor is selected if it is lower than the 
airspeed-oriented load factor. 


The pointing algorithm could generate negative. load 
factors. An option to command negative load factors in the 
maneuver-plane method has been partially implemented. The load 
factors commanded in the original AML were always positive. The 
equation for the maneuver plane is given by (p 53 of Reference 
6 ) : 


Rhos = arc tan 


-*e.T * x e,T *v 

(x e,T z e,T l x * ^e,T z e,T t y ‘ V h,T t z )/V T 


There are two solutions to this equation, Rhos and Rhos + 
180 degrees. The second solution corresponds precisely to 
negative a load factor, and is calculated in this version of 
AML. A negative load factor will be chosen if all these 
conditions are satisfied: 
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(1) B's airspeed must be lower than A's 

(2) A must be in B's forward quarter and low 

(3) B's current roll angle must not exceed 30 degrees; 
otherwise, it is preferrable to roll inverted under a positive 
load factor. 

(4) the negative load factor yields the smallest 
variation in maneuver plane rotation angle (and therefore in 
roll angle.) 

These conditions are restrictive and favor the well-known 
pilot preference for positive load factors. They will however 
make possible the use of a negative load factor for the purpose 
of bringing B's nose onto A while avoiding a high positive load 
factor and, hence, unnecessary loss of airspeed. 

Throttle Control The throttle control laws are set 
independently and can be summarized as follows: 

(1) In dive recovery, set the throttle to bring the 
airspeed near the corner velocity. Thus, the throttle is set to 
idle if the airspeed is above corner velocity. The throttle is 
otherwise set to afterburner. 

(2) Under other conditions, the avoidance of an 
overshoot takes precedence over the rule enunciated above. This 
will occur if A is in front of B and B has a high overtake 
velocity. In this case, the throttle is set to idle. 
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Figure 3-3. Summary of the AMLB Control Laws. 
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4. AIRCRAFT AND CONTROL SYSTEM DYNAMICS 


4.1 SELECTING AN APPROPRIATE MODEL 

What constitutes an appropriate model depends on the purpose of 
the simulation. As illustration, consider the two extreme cases: 
(1) Development of evasive maneuvers against an air-to air 
missile (2) Training of pilots in ECM tactics in a BVR 
environment. To capture the intricate dynamics between a highly 
agile missile and a fighter aircraft, it is necessary to 
simulate aircraft response to control surface deflections. This 
will rotate the aircraft in such a way that at any instant of 
time, the missile seeker head "sees" the aircraft under the 
proper aspect angle. In the BVR case, representing the aircraft 
as a point -mass may be adequate. Close-in visual air-to-air 
combat in a flight simulator lies somewhere between these two 
extremes. To achieve the necessary accuracy for the CIC 
simulation, two key performances of the aircraft must be modeled 
accurately: 

1) The Normal Acceleration 

2) The Roll Dynamics 

Normal acceleration determines how tight the fighter can turn 
and whether or not he loses energy during the turn. Roll 
performance determines how quickly the fighter can change the 
direction of the lift. In AML, roll performance is the 
determining factor in how fast the flight path can be changed 
from one maneuver plane to another maneuver plane. The two 
important parameters for roll performance are maximum roll rate 
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and maximum roll acceleration. As Shaw (Reference 10 ,page 414) 
points out: 

" in air combat, continuous rolls of more than 180 degrees 
are seldom required. Because a certain length of time is 
necessary to acclereate the roll rate from zero to its 
maximum value, maximum stabilized roll rate may not be 
reached during such short periods of roll. Therefore, roll 
acceleration is often the controlling factor in combat 
performance. " 


Shaw's quote is certainly true for air-to-air combat and even 
more so for evasive maneuvering against missiles. 


The problem of properly simulating roll performance is 
complicated by the fact that a change in bank angle has often to 
be achieved under high angle of attack or that coupled with a 
change in bank angle is a large change in angle of attack. In 
AML, a maximum roll rate and a maximum pitch rate is specified. 
Both are a function of the particular aircraft type represented 
by AML. If a maneuver command requires both a large change in 
the pitch angle (Theta hat) and the roll angle (Phi hat) the 
details of how this maneuver is performed depend a great deal on 
the ratio between maximum pitch rate and maximum roll rate. A 
proposed method, which, due to lack of funding never has been 
implemented, is to calculate the maximum available pitch and 
roll acceleration every time one of these extreme maneuvers has 
to be performed: 

As a first approximation. 
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is the 


rolling moment due to aileron deflection and do. 
maximum available aileron deflection. To be accurate, (' would 
have to be known as a function of Mach number and of the angle 
of attack. Herin lies the problem: It is often difficult to 

obtain these control derivatives for the extreme flight 
conditions which occur so often in air-to-air combat. Analogous 
remarks apply for C r. (control derivative for pitching moment 
due to elevator deflection) 

If the AML maneuver command is fed into a simulated (or 
eventually, into a real) flight control system, the problem of 
properly simulating pitch and roll performance under high angles 
of attack is greatly simplified. The aircraft (F-X) in 

Northrop' s Pilot Associate Program D1 was controlled by feeding 
AML provided load-factor and bank-angle commands into the flight 
control system. It can therefore be assumed that the dynamic 
response of the F-X to AML maneuver commands was very realistic. 

4.2 SYNOPSIS OF THE CURRENT ATTITUDE CONTROL MECHANISM 

A detailed account of the new equations of motion can be 
found in reference 14 where all the mathematical background 
underlying the treatement of the attitude control equations is 
presented. For* the sake of completeness of this report, the 
significant changes between the AML-75 and AML-84 are summarized 
below. 

As an introduction, a few words about "degrees of freedom" 
of an airplane model may be in order. If we consider the 
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aircraft to be a rigid body, then, by definition of classical 
mechanics, the number of degrees of freedom is equal to the 
number of independent coordinates required to uniquely define 
the position and the attitude of the body. A single rigid body 
can have at most six degrees of freedom (3 translational, 3 
rotational) . If we constrain the motion, the number of degrees 
of freedom is reduced, e.g. an aircraft whose eg could only move 
in a planeand whose longitudinal axis is constrqaint to lie in 
that plane, has 3 degrees of freedom (2 translational, one 
rotational) . How many degrees of freedom does the AML model 
have? The answer is this: We try to make it a five degree of 

freedom motion, by postulating that the sideslip angle and the 
rate of the sideslip angle (not the yaw angle and the yaw rate!) 
be zero. But during a transition from flight in one maneyuver 
plane into some other maneuver plane the calculated values of p 
q and r do not necessarily exactly guarantee a zero sideslip 
angle. The model is therefore a true six degree of freedom 
model . 

Most of the maneuver commands in AML are triplets defining 

- a maneuver plane ( by means of the maneuver plane 
rotation angle rho) 

- a load factor 

- a throttle setting 

Given the above three parameters, one can calculate what 
the aircraft's attitude, at the present time, (or one 
integration stepsize ahead) should be for the aircraft to fly in 
the commanded maneuver plane with the commanded load factor and 
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with zero sideslip. 


! 

Once we know the aircrafts desired attitude, we can 
calculate body rotational rates which will rotate the aircraft 
from its present attitude into its desired attitude. The 
important contribution of the "new equations of motion" is the 
way how these desired body rotational rates are calculated. To 
determine values of p q and r Euler angles Psi hat. Theta hat 
and Phi hat are calculated. These angles are expressed in the 
aircrafts present body axis system and not, as in the "old 
equations of motion" in the inertial reference system. 
Therefore, only Phi hat ever can become really large. Theta hat 
and Psi hat will always be relatively small ( Theta hat will 
never be greater than the difference between maximum and minimum 
allowable angle of attack) Consequently, there will never be a 
singularity in the set of Euler angles Psi hat. Theta hat and 
Phi hat, and as a consequence, the previously encontered problem 
of "going over the top " will no longer occur. 

4.3 REFINED CALCULATION OF COMMANDED PITCH RATE 

The procedure to determine p q and r as developed in 
reference 10 appeared to work reasonably well in the AML-84 
program, but occasionaly, the AML driven aircraft would fly into 
the ground even though dive recovery was inititiated at the 
appropriate time. Careful analysis of trajectories during dive 
recovery revealed that the aircraft never achieved the commanded 
load factor but consistently flew with a load factor less than 


4-5 


the commanded load factor during the pull out maneuver. At first 
we thought that the problem lies in the first order transfer 
function between q command and q achieved. But even as the time 
constant in this transfer function was reduced to a very small 
value, the problem persisted. The real reason for the 
discrepancy between commanded angle of attack and achieved angle 
of attack lies in the fact that, the calculation of the "desired” 
aircraft attitude is based on the present velocity vector. 
However, if the aircraft undergoes a large normal acceleration, 
the velocity vector will rotate during the next integration step 
and therefore, the commanded pitch rate must be increased by the 
rotational rate of the velocity vector which is: 



In a hard turn, a better value for q commanded therefore is: 

q com = min (q max , abs (0/dt)) sign (9) + Lift 

g*v 

4.4 THE POINTING CONTROL SYSTEM 

One of the most significant additions and improvements to 
the solution of the AML driven aircraft attitude control is the 
incorporation of a "pointing" control system. In several studies 
with AML, it was found that the AML controlled aircraft 
performed quite well to get behind the opponent, but once there, 
it lacked the capability to reduce the line of sight angle to 
the small value required for a gun solution. Controlling the 
aircraft by means of maneuver planes and loadfactors is indeed 
not a suitable way to point the aircraft's nose in a desired 
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direction. Therefore, a control system was implemented which 
would directly command roll and pitch rate to point the 
aircrafts longitudinal axis into a desired direction. Figure 4- 
1 illustrates in form of a block diagram the pointing control 
system. This control system is a modification of a control 
system suggested for use in surface-to-air missiles (Reference 
25, page 37). It is highly effective in controlling the AML 
driven aircraft. The problem is to find appropriate values for 
the various gains in the control system if a new fighter 
aircraft is implemented in AML. 
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p - FEEDBACK 


Figure 4-1 • Block Diagram of Pointing Control System 






5. SAMPLE TRIAL MANEUVERS AML VERSUS IF => THEN AML RUNS 


A series of test cases was conducted to exercise the 
AMLB logic described in section 3. In this series of runs, the 
'•A" aircraft was an F-15 controlled by the trial-maneuver logic 
AML (AMLA) . The "B" aircraft was an F-4, controlled by the IF => 
THEN AML logic (AMLB) . The initial conditions selected for these 
cases are shown in Table 5-1. A variety of initial velocities, 
altitudes and initial ranges were used. Initial velocities of 
M. 77 at 20,000 feet were selected because they represent a 
typical entry conditions into the air combat arena. On the 
other hand, initial velocities of M.46 correspond to typical 
corner velocities at 10,000 feet. The initial angular conditions 
vary from neutral to very unfavorable to the F-4 : the initial 
PI range approximately from 50 to 90. Since in addition the F-4 
is a considerably less performing aircraft than the F-15, which 
has a smaller turning radius, one would expect that the 
situation would develop in favor of the F-15. 
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Table 5-1 

Initial values for sample runs 


"A" —————— ii 3 h 


Region 

/Case 

Mach 

No. 

Altitude 

ft 

Relative 

Range 

ft 

Mach 

No. 

Altitude 

ft 

Eps(B) 

deg 

Lmbd(B) 

deg 

2/1 

.46 

10,000 

1000 

.46 

10,100 

90 

90 

3/1 

.46 

10,100 

1000 

.46 

10,100 

90 

135 

4/1 

.77 

20,000 

2250 

.77 

20,000 

154 

30 

4/2 

.77 

20,000 

2000 

.77 

20,000 

180 

45 

4/3 

.46 

10,000 

1000 

.46 

10,100 

135 

0 

4/4 

.46 

10,000 

1000 

.46 

10,100 

135 

45 

5/1 

.46 

10,000 

1000 

.46 

10,000 

135 

90 

5/2 

.46 

10,000 

1000 

.46 

10,000 

135 

90 

5/3 

.76 

16,000 

3500 

.77 

20,000 

135 

90 

5/4 

.77 

24,000 

3500 

.77 

20,000 

135 

90 

6/1 

.46 

10,000 

1000 

.46 

10,100 

135 

135 
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A natural way of classifying these runs is to define 
regions in the epsilon (B)- lambda (B) plane which is used for 
AMLB maneuver selection (Section 3) . This plane was divided in 
the regions shown below: 


LAMBDAB | 
degrees j 


180 


120 


60 


0 



EPSB, deg 


Figure 5-1. Definition of regions 
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These runs were made for a fixed period of time, typically 
20 seconds, which only allows the observation of the development 
of the initial maneuvers. The X-Y traces of the trajectories 
were plotted for that period of time. Also, additional pertinent 
information such as altitude and airspeed have been indicated as 
labels on these plots. For all these runs, the performance index 
for aircraft "A”, PI (A), and the integrated performance index, 
IPI(A), were plotted as a function of time. These were discussed 
in Section 1.6. Since PI(B) = 100 - PI (A), only PI (A) was 
plotted. The PI yields an indication of the relative angular 
attitude between the two aircraft and complements the 
information from the X-Y trace. It will be recalled that a PI of 
80 is required (but not sufficient) to achieve a firing 
position. 

5.1 REGION 2 

A run was made in this region corresponding to neutral 
conditions in all initial variables (angles, velocities and 
altitude.) These results are shown in Figures 5-2 and 5-3. 

It would be expected that the superior-turning F-15 would 
gain advantage. However, the performance index plot indicates 
that the situation remains essentially neutral throughout the 
simulated engagement. It will be noted that there is loss of 
airspeed on both sides (but is more severe for AMLA) , and that 
the engagement remains approximately co-altitude. 

An examination of the X-Y traces in Figure 5-2 contrasts 
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Figure 5-2 

Region 2 (A..(B)=90 deg., 
e(B)=90 deg.) X-Y Trace 
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4 5 6 7 8 9 10 11 12 13 14 15 16 


Time, sec 

□ PIB + IPIB 


Figure 5-3 Region 2 
PI Plot 




the initial behavior of AMLA/F-15 and AMLB/F-4E. AMLB begins to 
turn immediately, while AMLA reacts very little until t = 7 
seconds. Then AMLA begins a turn reversal during t = 7 seconds 
to 16 seconds. As indicated in Reference 10, this tactic should 
be expected from an aircraft with a smaller turn radius such as 
the F-15 compared to the F-4E. It will be seen again in more 
simulated runs. This tactic fails in the present case probably 
due to the small initial range (1000 ft) . 

The steering laws used by "B" have also been indicated in 
Figure 5-2. The pursuit law is used intermittently (t = 0 sec to 
1 sec; 12 sec to 13 sec; 14.5 sec to 16 sec) The lag pursuit law 
is used during the rest of the simulation. 

Both the "A” aircraft and the M B" aircraft rapidly lose 
airspeed, but remain approximately co-altitude during the 
simulated engagement. This trend will be observed in other 
engagements starting co-altitude and not involving tail-chase 
evasion. 
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5. 2. REGION 3 


The initial conditions used ( Lambda (B)= 135 degrees. 

Epsilon (B) = 90 degrees) put AMLB at a significant disadvantage 

with an initial PI(B) = 37. 

The trajectory, shown in Figure 5.4, exhibits flat scissors 
for the first 9 seconds of the engagement, during which the 
pursuit law is used. This part of the engagement is similar to 
the Region 2 case previously discussed. The F-15 has brief 
firing opportunities between t = 8 and 10 sec. Unlike the region 
2 case, AMLB cannot initiate a second scissor and is forced 
instead into the tail-chase evasion mode after t = 9 seconds, 

which accounts for its fluctuations in altitude. AMLB has a 
slight speed advantage during the major part of the simulated 
engagement. In a real-life engagement, this might be exploited 
to disengage, a maneuver not included in the present AMLB 
d i s engagement maneuver . 
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) X-Y Trace 





PERFORMANCE INDEX (A) 



TIME (SEC) 

□ PIA + IPIA 


Figure 5-5 Region 3 
PI Plot 
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5. 3. REGION 4 


Referring to Table 5-1, the range of initial epsilon(B)- 
lambda(B) values considered in region 4 corresponds to forward- 
quarter passes. This means that the initial position of each 
aircraft is in the other's forward quarter. As indicated in 
Reference 10 (p77) , there are two turn options available for 
fighters meeting in forward-quarter passes: the nose-to-nose 
turn option, and the nose-to-tail turn option. These are 
illustrated in Figure 5-6, adapted from Figure 2-11 in Reference 
10. The terminology refers to the position of the fighters at 
the end of the maneuver. These options were compared to the 
results obtained with the AMLA/AMLB logic in this series of four 
cases. 

Case 1 is illustrated in Figure 5-7 and shows a nose-to- 
tail conversion generated by the present AMLB. The instantaneous 
steering maneuver has also been indicated on Figure 5-7. "P" 
indicates that "B" follows a pure pursuit maneuver between t = 0 
seconds and t = 8 seconds. "LGP" indicates that a lag pursuit 
maneuver is used between t = 8 seconds and t = 21 seconds (end 
of the simulated engagement.) In this particular case, the AMLB 
steering law provides the F-4 with both a good defensive 
maneuver and a good maneuver for repositioning for attack. In 
contrast, the AMLA-controlled F-15 does not exhibit a 
repositioning tactic and seems instead to "wander off". The 
performance index plot (Figure 5-8) indicates that the situation 
evolves in favor of the F-4 from initially neutral conditions. 
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Figure 5-6 
Forward quarter 
turn options 











PERFORMANCE INDEX (A) 



Figure 5-8 Case 1, 
Region 4, PI Plot 
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Case 2 represents a slight variation in initial conditions 
compared to case is Lambda (B) = 45 deg; epsilon (B) = 180 deg. 
It is illustrated in Figure 5-9. In contrast to the previous 
run, this results in a nose-to-nose conversion. This conversion 
mode offers the F-4 the potential for a subsequent head-on 
firing opportunity when the range closes to less than 3000 feet 
(but this is beyond the interval simulated) . 

The maneuvers used in Case 2 have also been indicated on 
Figure 5-9. Pursuit (P) is steered between t = 0 second and t = 
1 second; between t = 2 seconds and t = 12 seconds; and between 
t = 14 seconds and 15 seconds. Lag pursuit (LGP) is steered 
between t = 1 second and t = 2 seconds; t = 12 seconds and 15 
seconds; and between t = 15 seconds to the end of the simulated 
engagement . 

The initial angular conditions for cases 3 and 4 also 
correspond to forward quarter passes, with slight variations in 
epsilon and lambda compared to cases 1 and 2. However, the 
initial altitudes (10,000 feet) and speed (M.46) are very 
different. The initial speed was selected so that both aircraft 
start near corner velocity, the velocity at which both aircraft 
have their best turn performance. The X-Y traces for both cases 
(Figures 5-11 and 5-13) rapidly develop into well-defined 
“scissors". The effect of the F-15's smaller turn radius is 
apparent: "A" turns well within "B". However, in spite of this 
visible advantage, the PI plots for both cases indicate that "A" 
does not attain a gun-firing position. The situation remains 
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Figure 5-10 
Case 2, Region 4, 
PI Plot 
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essentially stalemated. 


The steering laws have also been indicated on Figure 5-11 
and 5-13. The pursuit law is used the most frequently, as would 
indeed be expected from the domains specified in Figure 3-2. 


3000 













Figure 5-12 
Case 3 PI Plot 
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Figure 5-14 
Case 4 PI Plot 


5-22 



5.4. REGION 5 


A series of 4 runs were made in this region which 
corresponds to a slight angular advantage in favor of the F-15. 
In case 1 of region 5, shown in Figure 5-15, the F-15 eventually 
gains angular advantage after t = 13 seconds, but also loses 
more airspeed than the F-4E in the turn. 

Case 2 of region 5 was run to highlight the influence of 
the turning ability of the F-15 on the result of an engagement 
with the same initial conditions as in Case 1. The thrust/weight 
ratio of the F-15 was reduced by increasing the weight from the 
nominal 40,000 lbs to an artificial 50,000 lbs, thus yielding a 
thrust-to-weight ratio of approximately 0.8 which is comparable 
to the F-4E . The results are illustrated in Figures 5-17 and - 
18. This run shows that the F-4E now has a firing opportunity 
between t =11 and 12 sec. 

Cases 3 and 4 illustrate the effect of an initial altitude 
difference on the same initial angular conditions as in Case 1. 
In case 3, "A M has an initial altitude advantage of 4,000 feet 
compared to "B". The results have been illustrated in Figures 
5-19 and -20. The present AMLB logic commands a pursuit course 
with the aimpoint located at the altitude of "A”. As a result, 
"B” rapidly loses altitude. The situation at the end of the 
simulated engagement shows that "B" ends up in a defensive 
position. Thus, the initial altitude advantage has not improved 
"B"'s situation. 
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Figure 5-13 Region 5, 
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Figure 5-18 Region 5, 
Case 2, PI Plot 
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The PI at the end of the simulated engagement ( Figure 5- 
20) was much higher than one would expect at first from an 
examination of the X-Y trace. For this reason the line-of-sight 
and angle-off were plotted individually in Figure 5-21. A 
careful examination of the run shows that the "A” aircraft is 
strongly pitched down. This attitude explains the observed 
variation in these angles. This highlights the utility of the PI 
in summarizing the angular situation of the engagement. 

Case 4 assumes an initial altitude disadvantage of 4,000 
feet for "B". In this case, illustrated in Figure 5-22, the AMLB 
logic commands a climbing turn in "A'"s direction, resulting 
from the pursuit law which is used between t - 0 seconds to 7.5 
seconds. This maneuver brings "B" in "A” 's forward quarter, but 
the PI plot in Figure 5-23 shows that "A” does not have a firing 
opportunity as a result of the altitude difference. The 
engagement ends up with "A" overshooting "B", without "B" having 
a gun-firing opportunity due to the altitude difference. 
Following the overshoot, the F-15 does not appear to be 
reacting. 

In both cases 3 and 4, both "A” 's and "B" 's tactics could 
be improved by the inclusion of negative G's. 
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Figure 5-21 Region 5, Case 3 
Line-of-Sight and Angle-off Plot 
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5.5. REGION 6 


This region corresponds to a severe initial angular 
disadvantage for "B". The case illustrated in Figures 5-24 and 
5-25 exercizes the AMLB evasive maneuver during the entire 
simulated engagement. The PI plot indicates that there is no 
improvement in "B' M s angular position. However, the X-Y trace 
indicates that the relative range increases from an initial 1000 
feet to 3500 feet. This is due to "A" ' s rapid loss of airspeed 
during the turn. This result shows suggests that "B H might have 
an opportunity to disengage. 

5.6. CONCLUSIONS 

In all cases, a wide difference in outcomes has been 
observed for small variations in the initial angular conditions. 
This result has often been observed in ACM simulations. 
Furthermore, this wide difference in outcomes occurred in spite 
of a small variation in initial PI. The use of the initial PI to 
classify and predict the entire engagement outcome does not 
appear promising. 

In all cases, both AMLA and AMLB command initially high 
load factors which result in a rapid loss of airspeed. In all 
cases starting co-altitude, the fight remains roughly in the 
initial horizontal plane unless tail-chase evasion is initiated. 
The lag-pursuit and pure pursuit laws involve essentially a 
series of level turns, or an "angles fight" to use the 
terminology of Reference 10. In this fight, the F-4 cannot gain 
an advantage due to its lower turning capability compared to the 
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Figure 5-25 Region 6 
PI Plot 
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F-15. In many cases, surprisingly, the fight does not rapidly 
evolve to the disadvantage of the F-4, as might be expected from 
the disparity in performance, and remains approximately neutral. 
Finally, it will be noted that the cases investigated did not 
present any opportunity to exercize the "pointing algorithm" 
dsicussed in Section 3. 

Throughout the history of air combat, skilled pilots have 
been able to win engagements in spite of having the lower- 
performing aircraft. In this situation, they would avoid a 
turning fight as simulated above. In the next section, we 
describe an alternative approach to angular conversion which 
attempts to trade off altitude to gain an angular advantage. 
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6. ADDING BASIC FIGHTER MANEUVERS TO THE IF => THEN LOGIC 


The maneuvers generated by the IF => THEN AMLB and 
discussed in Section 5 are realistic and have generally enjoyed 
good pilot acceptance in flight simulators. However, the 
predictability of the maneuvers it generates has been criticized 
because it makes it possible for a person to anticipate AMLB ' s 
future maneuvers after a few sessions in a simulator. To enrich 
the variety of maneuvers generated by AMLB, additional maneuvers 
based on the "basic fighter maneuvers'* (BFMs) of the type found 
in ACM training manuals (for example, references 10, 11, 12) were 
added to the existing AMLB logic. It will be recalled from 
Section l that such an approach had been rejected at the time of 
the development of the original AML program. However, BFMs were 
used in the present effort because they improve the variety of 
maneuvers generated by AMLB, not only in flight simulators, but 
also against AMLA in offline programs, and has proven useful in 
these respects. 

In examining samples of such BFMs, it was found that in 
general each of these maneuvers is appropriate under a narrow set 
of circumstances based primarily on relative geometry, and 
additionally on other situational parameters such as closing 
velocity, relative airspeed, relative altitude, to name just a 
few. While all these maneuvers have their individual, specific 
objective, the majority of them share with AMLB the underlying 
purpose of angular conversion on the opponent, except in the case 
of the disengagement maneuver, which does not exist in AMLB. This 
observation suggested keeping the underlying AMLB angular 
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conversion logic, and to replace it only when an opportunity for 
a BFM arises. To reduce predictability, when a choice between a 
BFM and AMLB is available, the selection of the maneuver is 
decided by the means of a random number. 

Due to the limited scope of this work, only a set of three 
BFMs was investigated: a "diving overshoot", a "vertical 
overshoot", and "opposite turn". This terminology and the results 
obtained will be discussed in detail later in this section. 

The features discussed above were implemented in a new 
subroutine called SELECTB. The particular requirements for the 
BFMs were (1) to identify when a particular BFM can be executed; 
(2) to execute that maneuver for a specified amount of time, or 
until conditions specific to that BFM are no longer met, and (3) 
to terminate the maneuver under the specified conditions and 
return control to the underlying AMLB. Function (1) is presently 
performed in a added subroutine called SELECTB, while functions 
(2) and (3) are performed in individual subroutines, as 
illustrated in Figure 6-1. 
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SUBROUTINE SELECTS 


| Ground Avoidance Test | 

+ — + 

I +— 

I — YES — | Execute Ground Avoidance I > RETURN 

NO | 


+ + 

I BFM in Process? | 



I 

| YES 

NO 


+ ------- — --- — --- t- 

| Identify all possible | 

I BFMs based on relative j 

j geometry, speed, etc.. | 

I I 




■+ 


| BFM subroutine | 
— > j library j 

+ + 


I 

| - — BFM in process > RETURN 

No BFM in process j 


v v 

| Execute BFM selected | 
I Call appropriate subroutine j 
j from BFM library j 



Figure 6-1 OUTLINE OF THE IMPLEMENTATION 

OF BASIC FIGHTER MANEUVERS 


6-3 




No BFM in process | 



j Pointing algorithm decision | 



V 

REACTB logic 


+ + 


j Lead/Lag algorithm decision 




j Load Factor determination I 
| (Positive/Negative) j 





I 

j Evasive maneuver decision 

I 





RETURN 

Figure 6-1 Concluded 
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6.1 OPPOSITE TURN 


Several cases of "forward passes" were discussed in Section 
5. It was found that the AMLB logic would generate either a nose- 
to-nose conversion, or a nose-to-tail conversion, depending on 
the initial conditions. Reference 10 (p 79) outlines the 
potential advantage of a nose-to-nose turn, which could result in 
achieving offensive advantage. To execute such a maneuver in Case 
1 of section 5 would require "B" to turn away from his opponent, 
as was illustrated in Figure 5-1. This maneuver was 
implemented in SELECTB as the "opposite turn" BFM. 

Sample results are shown in Figure 6-2 and 6-3. The effects 
of the superior turning ability of the F-15 over the F-4 is 
clearly demonstrated by the AMLA-controlled F-15 performing a 
"turn reversal" to gain a firing position on the F-4. This 
possibility of this situation is in fact predicted in Reference 
8, Figure 2-12. In this case, AMLA derives the "textbook 
solution". The F-4 does not gain anything by performing an 
opposite turn with the initial conditions considered. 

To underscore the role of the turn rate, another case was 
run with the same initial conditions as in case 1, but now with 
the F-E initially at corner velocity (M .41), and the F-15 
remaining at the same initial velocity of M .77 . Case 2 has been 
illustrated in Figures 6-4 and 6-5. The F-E clearly achieves a 
nose-to-nose offensive position against the F-15. The F-15 also 
achieves an offensive position. 

The results of these cases would encourage us to consider 
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the differential turn rate between the two aircraft at the time 
of the decision for the next maneuver. The implementation of 


this feature would be fairly simple from a 
standpoint. Furthermore, AML presently does not 
between a "defensive" aircraft and an "offensive" 
inclusion of the differential turn rate might be a 
so. 
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aircraft. The 
good way to do 
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Figure 6-5 Opposite turn. 
Case 2 PI Plot 
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6.2 DIVING OVERSHOOT 


This BFM is a gun defense maneuver described in Reference 10 
(pages 26-27) . The defending aircraft dives in order to force the 
opponent to overshoot. For brevity, it is referred to in this 
document as a "diving overshoot". It is activated under the 
following conditions: 

1. 120 deg <■ Epsilon (B) <= 180 deg AND 

2. 60 deg <« Lambda (B) <= 120 deg AND 

3. ABS(ZEA - ZEB) <■ 1000 feet (i.e. approximately 

co-altitude) , AND 

4. Relative Range <* 6000 feet 

It is implemented in subroutine OVRSHT. The diving effect is 
accomplished by steering the B aircraft on a pursuit course which 
uses an aimpoint 10,000 feet below the A aircraft, i.e. with 
coordinates XEA, YEA, ZEA + 10000 (z positive downwards) . The 
commanded load factor is 95% maximum. In spite of the dive, the 
commanded throttle setting is A/B because of the anticipated 
speed loss due to the high-G turn. Due to the anticipated loss of 
altitude, the maneuver can ony be executed above a minimum 
altitude. 

This maneuver was only tested against a non-interacting "A" 
aircraft, i.e. flying straight and level. The results are shown 
in Figures 6-6’ and 6-7 (Case 1). A case with the original AMLB 
lead/lag logic is illustrated in Figures 6-8 and 6-9 (Case 2) . 
The X-Y traces are dramatically different. However, there is less 
difference than anticipated in the PI plots. In Case 1, it will 
be observed that the "diving overshoot" results in a 5700 foot 
altitude drop for "B" . "B" 1 s velocity decreases in spite of the 
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altitude drop and full afterburner setting. It can also be 
observed from Figure 6-6 that a large portion of the altitude 
drop occurs after the "diving overshoot command" (which is 
indicated on the figure as "DOVS") is replaced by the 
conventional pursuit (indicated as "P") after t * 3 seconds. 

The maneuver is terminated when the range rate increases. 
However, an examination of case 1 suggests that the maneuver 
might instead be terminated earlier to avoid the altitude drop. 
For example, a criterion for maneuver termination might be when 
"B" has crossed "A"'s projected track. 
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O n A" (F-15) Figure 6-8. Lead/Lag overshoot 

(F-4) X-Y trace 
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Mach No. /Altitude (in boxes) 



PERFORMANCE INDEX (A) 


100 
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6.3 PULL-UP OVERSHOOT 


In section 5, some conditions under which scissors maneuvers 
are generated in forward passes were described. Reference 19 
(page 6-28) describes a counter to that maneuver which consists 
initially of a pull-up with the intent of an overshoot. In order 
for this maneuver to work, the "B" aircraft needs a velocity 
excess which it can convert into an altitude advantage. This 
maneuver will be referred to here as a "pull-up overshoot". 

A sample case is shown in Figure 6-10 and 6-11. These are 
neutral initial angular conditions, with the "B" aircraft having 
a speed advantage over the "A" aircraft (M.74 vs M.46). The 
initial conditions are similar to Case 4 of Section 5 shown in 
Figure 5-8. 

The maneuver consists of a wings-level high-G (95% of 
maximum G) pull-up for a specified number of seconds (in this 
case 5 seconds.) Following the pull-up, control is reverted to 
the lead-lag logic to finalize the angular conversion with more 
favorable parameters. The wings-level pull-up causes the F-4 to 
separate angularly from the F-15. Because of the relative 
position at t- 5 sec (the time the F-4 stops climbing) , the next 
maneuver was a tail-chase evasion (indicated in the figure as 
"TCE" ) , followed finally by a pursuit ( indicated as "P") . As 
shown in Figure 6-11, the trend indicated by the performance 
index is that "B" is gaining an angular advantage. This would 
indicate that the maneuver has proven succesful in breaking the 
stalemate of the scissors in "B" ' s advantage. 








7. CONCLUSIONS AND RECOMMENDATIONS 


In the almost twenty years since work on the adaptive 
maneuvering logic started, the AML programs have been 
continuously, little by little, improved. The AML version, which 
was distributed around 1978 by COSMIC had severe deficiencies, 
mostly in the attitude dynamics. It was unfortunately this 
versions (or derivatives thereof) which were installed in a 
number of US Navy and Air Force flight simulators. 

Todays state of AML is that the motion of the AML driven 
aircraft is quite realistic and AML's tactical behavior is most 
of the time sound. 

In the course of this long development period, we have 
learned a few basic and important lessons for the simulation of 
close-in air-to-air combat and for missile evasion: 

- realistic aircraft motion, specifically the 
rotational dynamics, is of greatest importance for 
pilot acceptance. 

- accurate roll and pitch dynamics are crucial when 
developing evasive maneuvers against surface-to-air 
or air-to-air missiles. 

- improving the tactical behavior of AML is very time- 
consuming and tedious. 

- the performance of any air-to-air combat program can 
only be evaluated statistically. Well over 100 
different initial conditions must be exercised to 
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arrive at valid statistics. 

both methodologies, trial maneuver and IF => THEN, 
show promise. 

- an analysis of the performance of an air-to-air 
combat program can not be made by analysis of non- 
real time, batch processing runs alone. An 
interaction with highly skilled human pilots is 
absolutely required. 

- real-life air-to-air combat is extremly complex. 

The original idea of this contract was to prepare a real- 
time base-line version of AML which could be used by the Flight 
Research Center to play the role as a "flight-director” 
controlling an actual airborne aircraft. By uplink telemetry, 
maneuver commands are issued to the aircraft and by down-link 
telemetrty, aircraft status is received. Thus, the entire 
computational effort can be performed on ground. The complexity 
of such a project precluded implementation under this contract. 
We did, however, provide the Flight Research Center with an IF 
=> THEN version of AML, running in real time in conjunction with 
an existing flight simulation. Due to lack of adequate real time 
display facilities, this air-combat simulation was not used much 
by the Flight Research Center. 

We also recognize, at this point, that a number of problems 
in the simulation of one-versus-one combat still require 
additional studies and analyses. To name just a few: 
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- How can we prevent, early enough, the AML driven 
aircraft's energy to deteriorate to a very low value? 

- If the AML driven aircraft "flies” against a dissimilar 
aircraft, how do we make best use in performance 
differences between the two aircraft? 

- Is it possible to build an AML where all the IF => THEN 
production rules are formulated in plain English, so 
that a fighter pilot can change them at his will and 
investigate the effects of the change? 

- How can distributed and parallel processing help to 

overcome some of the limitations imposed presently on 
real-time versions of AML? 

- How can we put AML on-board a remotely piloted aircraft 

and then perform the ultimate "proof" for AML's tactics? 
How can we incorporate some of the results of 
the theory of differential games into AML? 

If one admits that the decision logic of AML is not yet 
perfect (and the authors of this report certainly admit that) , 
then a challenging problem is the following: How can we 

methodically improve AML? There are two aspects: (1) How to make 
changes to the decision logic and (2) How to evaluate the 

effects of these changes. It appears that a solution to this 

problem requires extensive use of a real-time full-dome flight 
simulator and the cooperation of experienced fighter pilots. 

It appears, that after almost twenty years, the challenges 
in building an "iron-pilot" have not become smaller, but have 
grown. 
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APPENDIX 

LISTING OF THE FORTRAN ROUTINES FOR 
"BASIC FIGHTER MANEUVERS" 


SUBROUTINE 

SELCTB 

A-l 

SUBROUTINE 

OVRSHT 

A-8 

SUBROUTINE 

VTOSH 

A-10 

SUBROUTINE 

OPSTRN 

A-12 


SUBROUTINE CLIMB 


A- 14 


1 7— Ju 1 - 1987 14:0’ 
1 7- Ju 1-1 987 14:0: 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
001 1 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 


SUBROUTINE SELCTB (XEA, YEA* ZEA, XEDOTA, YEDOTA, ZEDOTA, DMTRXA ) 

C 

C 

C SUBROUTINE BASED ON REACTB (MAR 1986) TO SELECT AND EXECUTE 

C APPROPRIATE MANEUVERS. WHEN SEVERAL MANEUVERS ARE FEASIBLE* ONE 

C SPECIFIC MANEUVER IS SELECTED USING A RANDOM NUMBER 

C WHEN NO SUCH MANEUVER IS FOUND, THIS SUBROUTINE REVERTS TO THE 

C LEAD/ LAG STEERING LAW OF REACTB 

C 

c 

C 

C *** COMMON BLOCKS FOR REAL TIME SIMULATION AT NASA DRYDEN 
C 

COMMON/ CNSTNS/ DT , TBEGN, TNOW, P I , P I DV2, PI DV4, TWOP I , DEGRD , RADDG, G, 

1 VAR (20) , I VAR (20) , TEND 

C 

COMMON/ CONTRL/MSTOP, I PR I NT 
C 

COMMON/ COMNDB/ I CMNWB , GLEVLB , ROTB , MANVRB 
C 

COMMON/ RNDM AN/ MLDECS, MAN I N I , I SLCTR 

COMMON/ DATA 1 B/ XEB, YEB, ZEB, XEDOTB, YEDOTB, ZEDOTB, XEDDTB , YEDDTB , 

1 ZEDDTB, PSIB, THETAB, PH IB, UB, VB, WB, PB. QB, RB, A1B, A2B , 

2 A3B , A4B, VELB , VHORB 
C 

COMMON/ DAT A2B/ ALFAB , BETAB , CBARB (3,3) , CDB , CLALFB , DMTRXB (3.3), DRAGB , 

1 L I FTB , LODMXB , LODSTB , MACHB , RHOB , SPECEB , SB , THRSTB , 

2 PSUBSB , TPOSB » I N I ZB . WE I TB , CSB , CLB , PS I BRB , THETBB , 

3 AN1B, AN2B , AN3B , MASSB 
C 

REAL L I FTB , LODMXB , LODSTB , MACHB , MASSB 

C 
C 
C 

COMMON/ TBF4EB/CLMAXB ( 14) ,XTAB1 (14) ,NX1, 

THR I DB (7, 14) , XTAB2 ( 7 ) ,YTAB2(14) , NX2 , NY 2 , 

THRMLB <7, 14) ,XTAB3(7) ,YTAB3(14) .NX3.NY3, 

THRABB (7, 14) ,XTAB4(7) ,YTAB4(14) ,NX4,NY4, 

ALFCLB ( 16,10) ,XTAB5( 16) , YTAB5 ( 10) , NX5 , NY5 , 

CLFALB( 10, 11), XTAB6( 10) , YTAB6 ( 1 1 ) ,NX6, NY6, 

CDFCLB <16, 10) , XTAB7 (16) ,YTAB7(10) ,NX7,NY7, 

CLFCDB (18, 10) , XTAB0 (18) ,YTAB8(10) ,NX8,NY8, 

RECAGB < 10 , 12) , XTAB9 (10) , YTAB9(12) ,NX9,NY9 

C 

C 

COMMON/ PO I NTP/ XEA I M , YEA I M , ZEA I M , PCOMB , QCOMB , RCOMB , I PO I NT 
C 

DIMENSION CMPL (3,3) 

C 

COMMON/RELVAR/LOSELA, LOSELB, LOSAZA, LOSAZB, LOSANA, LOSANB , 

1 LSDOTA, LSDOTB, DEVANA, DEVANB, DVDOTA, DVDOTB, 

2 RANGE , RRATE , XA I NB , YA I NB , Z A I NB , XB I NA , YB I N A , 

3 ZBINA, ANGOFA, ANGOFB 
C 

REAL LOSELA , LOSELB , LOSAZA , LOSAZB . LOSANA , LOSANB , 


A-l 


1 


LSDOTA, LSDOTB 


SELCTM 


0058 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 
0081 
0082 

0083 

0084 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 
0093 

0099 

0100 
0101 
0102 

0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 
011 1 
0112 

0113 

0114 


1 7— Ju 1 — 1987 
1 7- Ju 1-1 987 


14:0 

14:0 


C 

C 

C *** LET THROTTLE ROUTINE KNOW WHEN A/C IS IN DIVE RECOVERY 
C 

COMMON/ D I VEB/ I RECVB 
C 

COMMON/PAGECT/ I CNT 
C 

DATA I XL 1 » J YL 1 * IXL2, JYL2/ 1 , 1,1,1/ 

DATA IXL9, JYL9/1, 1/ 

C 

C THE DEFAULT DECISION INTERVAL IS MLDEF , ELSE IT IS SET IN 

C THE INDIVIDUAL MANEUVER ROUTINE 

C 

DATA MLDEF/ 20/ , NEGTVG/ 1 / 

C 

IF (INIZB.EQ.l) THEN 
C 

ML DECS= MLDEF 
I SLCTR=0 
I RECVB=0 
RETURN 

END IF 
C 


C GROUND AVOIDENCE LOGIC 

C 

c 

c 

HB=-ZE8 

D I VEAN=-THETBB 

I F ( I RECVB . EQ . 1 .AND. D I VEAN . GT. 0. ) THEN 
MANVRB= 1 
I CMNWB= 1 
ROTB=0 . 

I PO I NT =0 
GO TO 998 
END IF 

IF< IRECVB. EQ. 1 .AND. D I VEAN . LE. 0. ) THEN 
I RECVB=0 
ELSE 

J F ( HB. LT. 20000. ) THEN 
FMACHX=MACHB 

I F ( FMACHX. LT. . 4) FMACHX=. 4 
HX=HB 

I F ( HX. LT. 200. >HX=200. 

CALL TLU2<HX, FMACHX, XTAB9, YTAB9, RECAGB , NX9 , NY9 , 
I XL9, JYL9, RECAN, IC) 

ELSE 

RECAN=P I DV2 
END I F 

I F ( D I VEAN . GT . RECAN ) THEN 
ROTB=0 . 

GLEVLB= 1 . 


A- 2 



SELCTll 


17-JU1-1987 14:0 
1 7- Ju 1 - 1987 14:0 


0115 

0116 

0117 

0118 

0119 

0120 
0121 
0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

0157 

0158 

0159 

0160 
0161 
0162 

0163 

0164 

0165 

0166 

0167 

0168 

0169 

0170 

0171 


C 

C 

C 

c 

c 

c 

c 

c 

c 


c 


c 


c 

c 

c 

c 

c 

c 


c 


c 

c 

c 

c 

c 


ICMNWB=1 
I PO I NT-0 
IRECVB*! 
MANVRB«2 
GO TO 998 

END IF 
END IF 


*** END OF GROUND-A VO I DENCE LOGIC 


2222222 START OF SELECTABLE MANEUVERS 22222222222222222222222222: 


IF ( MOD ( I VAR ( 1 ) * MLDEF ) . EQ. 0 .AND. ISLCTR .EQ. 0) THEN 
EPSD=ANGOFB*DEGRD 
FLAMBD*LOSANB#DEGRD 
MAN I N I =— 1 



IF 

(EPSD .GE. 

120. 

.AND. EPSD .LE. 

180. 

. AND. 

1 


FLAMBD. LE. 

60. 

) THEN 







I SLCTR- 101 

0 



END IF 






IF 

(EPSD .GE. 

120. 

.AND. EPSD .LE. 

180. 

. AND. 

1 


FLAMBD. GE. 

60. 

.AND. FLAMBD. LE. 

120. 

) THEN 


IF ( ABS < ZEA-ZEB ) .LE. 1000.0) THEN 


OPTION 1 (DEFENSIVE) FORCE AN OVERSHOOT 
OPTION 2 (DEFENSIVE-OFFENSIVE) ROLL REVERSAL 

I F (RANGE .LE. 6000.) ISLCTR=400 

I SLCTR-200 


ELSE IF (ZEA . LT. ZEB ) THEN 
ISLCTR=810 

ELSE IF (ZEA .GT. ZEB) THEN 

ISLCTR*910 

END IF 

END IF 


IF (EPSD .GE. 120. .AND. EPSD .LE. 180. .AND. 
1 FLAMBD. LE. 60. ) THEN 

1 SLCTR»6 1 0 
END IF 

WR I TE ( 77 » 491 ) TNOV , I SLCTR * EPSD . FLAMBD , RANGE 

WRI TE ( * * 491 ) TNOW , I SLCTR * EPSD* FLAMBD > RANGE 

% 


END IF 
C 
C 


ISLCTR*0 

C 

C 
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0172 

C 

2222222 EXECUTION PART OF SELECTABLE MANEUVERS 

0173 

0174 

c 

IF 

(ISLCTR - NE. 0 .AND. MOD ( I VAR ( 1 ) , MLDFF ) . EQ . 0 ) THEN 

0175 



I F ( ISLCTR . EQ. 400) 

0176 

0177 

0178 

c 

1 

CALL OVRSHT ( XEA, YEA , ZFA * XEDOTA, YEDOTA, ZEDOTA, DMTRXA ) 

0179 



I F( ISLCTR .EQ. 810) 

0180 

0181 

c 

1 

CALL VTOVSH ( XEA , YEA * 2EA , XEDOTA, YEDOTA, ZEDOTAf DMTRXA ) 

0182 

0183 

c 


I F ( ISLCTR . EQ. 610) 

0184 

0185 

c 

1 

CALL. OPSTRN ( XEA. YEA* ZFA. XEDOTA, YEDOTA, ZEDOTA, DMTRXA) 

0186 



I F ( ISLCTR .EQ. 1010) 

0187 

0188 
0189 

c 

1 

CALL. CLIMB (XEA, YEA, ZEA, XEDOTA , YEDOTA, ZEDOTA , DMTRXA ) 

0190 


END IF 

0191 

c 



0192 

0193 

c 

IF 

(ISLCTR -NE. 0 ) RETURN 

0194 

c 

2222222 END OF SELECTABLE MANEUVERS 222222222222222222222222222: 

0195 

c 



0196 

c 

START OF LEAD/ LAG 111111111111111111111111111111111111111111111111. 

0197 

c 



0198 

c 

MLDECS=200 

0199 


IF 

( MOD ( I VAR ( 1 ) , MLDECS ) . EQ . 0 ) THEN 

0200 



I PO I NT=0 

0201 



EPSD=ANGOFB*DEGRD 

0202 



FLAMBD=LOSANB#DEGRD 

0203 

0204 

c 


ROTPRV-ROTB 

0205 

c 

* ## 

CHECK FIRST IF WE WANT TO INVOKE THE POINTING ALGORITHM 

0206 

0207 

c 


I F ( FLAMBD. LE. 30. .AND. EPSD.LE.45.) THEN 

0208 



I PO I NT= 1 

0209 



XEA I M=XEA 

0210 



YEA I M=YEA 

0211 



ZEA I M=“ZEA 

0212 



MANVRB=7 

0213 



GO TO 998 

0214 

0215 

c 


END IF 

0216 

c 

# ## 

LEAD-LAG -PURSUIT DECISION FOLLOWS 

0217 

0218 

c 


IF (EPSD. LE. 30. .AND. FLAMBD. LE. 30. ) THEN 

0219 



DTPRED-3. 

0220 



MANVRB=3 

0221 



ELSE IF (FLAMBD. LE. (90. -EPSD) ) THEN 

0222 



DTPRED=0. 

0223 



MANVRB=4 

0224 



ELSE I F ( FLAMBD. LE. (180. -EPSD)) THEN 

0225 



DTPRED=-3. 

0226 



MANVRB-5 

0227 



ELSE 

0228 



DTPRED=0. 
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0229 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0238 

0239 

0240 

0241 

0242 

0243 

0244 

0245 

0246 

0247 

0248 

0249 

0250 

0251 

0252 

0253 

0254 

0255 

0256 

0257 

0258 

0259 

0260 
0261 
0262 

0263 

0264 

0265 

0266 

0267 

0268 

0269 

0270 

0271 

0272 

0273 

0274 

0275 

0276 

0277 

0278 

0279 

0280 
0281 
0282 

0283 

0284 

0285 


MANVRB=9 
EN1) IF 

XEXA=XFA+DTPRED#XEDOTA 
YEXA=YEA+DTPRED# YEDOTA 
ZEXA=ZEA+DTPRED#ZEDOTA 

TAXE=XEXA-XEB 
TAYE=YFXA-YEB 
TA7.E=ZEXA-ZEB 
VH0R2=XED0TB##2+YED0TB* #2 
VHORB=SQRT ( VH0R2 ) 
VEL.2=VH0R2+ZED0TB**2 
VEI .B=SQRT ( VEL2 ) 


C 

C 

C 


c 

c 

c 

c 


c 

c 

c 

c 


DZ" ( XEDOTB»ZEDOTB*TAXE+YFDOTB*ZEDOTB#TAYE-VHORB* *2#TAZE) / VELH 
DY- -YEDOTB*TAXE+XEDOTB*TAYE 
I F ( DZ. EQ . 0. .AND. DY.EQ.0.) THEN 
ROTB=0. 

ELSIE IF (DY.EQ.0.) THEN 
IFCDZ.GT.0.) ROTB=0. 

IFCDZ.LT. 0. ) ROTB=P I 
ELSIE 

R0TB=ATAN2 ( DY * DZ ) 

END I F 

*** SELECT THE POSITIVE G-LEVEL DEPENDING ON B’S VELOCITY 

I F ( VELR.GT.400. ) THEN 

GLVPOS= ( LODSTB+LODttXB ) / (2. #LODMXB ) 

ELSIE 

GLVPOS=LODSTB/ LODMXB 
END IF 

T AS = VELB#0. 5925 

CAS=TAS #SQRT ( RHOB/ 0 . 0023768 ) 

IFCCAS.GT. 1.2*330. ) GLVPOS= ( LODSTB -1 LODMXB ) / (2. #LODMXB ) 

I CMNWB=1 

*** CALCULATE INTERCEPT TRAJECTORY G-LEVEL 

I F 1 ( FLAMBD. LT. 60. ) THEN 

CALI. D I RCOS ( PS I BRB * THETBB , ROTB , CHPL ) 

D I ST2-TAXE# #2+TAYE* *2 + TAZE* *2 

ZMT-TAXEfcCMPL (3# 1 ) +TAYE*CNPL (3* 2 H TAZE*CMPL < 3, 3) 

RAD I S~D I ST2/ ( 2 • #ZMT) 

GL2=* ( ABS ( ( VELB# #2 ) / RAD I S ) / G ) + CMPL (3,3) 

GL3~ABS (CMPL (2,3) ) 

GLEVRB=SQRT<GL2*#2+GL3* *2) /LODMXB 
I F ( GLEVRB . LT. GLVPOS ) THEN 
GLVPOS=GLEVRB 
END IF 


CALCULATE INTERCEPT TRAJECTORY FOR NEGATIVE G’S 


21 APRIL 


I F < ROTB .LE. 0. ) THEN 
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0287 
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0289 
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0291 
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0305 

0306 

0307 

0308 

0309 

0310 
031 1 

0312 

0313 

0314 

0315 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


ROTBNG-PI+ROTB 

ELSE 

ROTBNG=-P I +ROTB 

END IF 

CALL D I RCOS ( PS I BRB * THETBB * ROTBNG . CMPL ) 

ZMT=T AXE* CMPL <3# 1 ) +TA YE*CMPL ( 3, 2 H TAZE*CMPL (3*3) 

RAD I S- D I ST2/ (2. *ZMT) 

GL2= ( ABS ( ( VELB# *2 ) / RAD I S ) / G ) +CMPL (3,3) 

GL3—ABS ( CMPL (2,3) ) 

GLEVRB=SQRT (GL2**2+GL3* *2 ) /LODMXB 
GLVNEG=2. 0/ LODMXB 
I F ( GLFVRB . LT. 2 . 0/ LODMXB ) THEN 
GLVNEG=GLEVHB 
END IF 

DETERMINE IF NEGATIVE G’S ARE ALLOWED, INCLUDING WHETHER GLVNEG EXCE 
ALLOWABLE LEVEL (PRESENTLY SET TO -2 G) 

CONDITIONS FOR US I NG NEGATIVE G f S: 

1. B’S AIRSPEED MUST BE LOWER THAN A’S 

2- A MUST BE IN B’S FORWARD QUARTER AND LOW ( FLAMBD <60 DEG, 

LOSELA < -5 DEG) 

3. B’S ROLL ANGLE MUST NOT EXCEED 30 DEGREES- ELSE HE IS BETTER OFF 
ROLLING INVERTED 

THE OBJECTIVE OF NEGATIVE G’S IS TO BRING B’S NOSE ON A WITHOUT 
PULLING HIGH G’S (HENCE LOSING AIRSPEED) 

(A FUNCTION HAVING SIMILARITIES WITH THE POINTING ALGORITHM) 

IN GENERAL., NEGATIVE G f S WOULD BE USED TO UNLOAD THE AIRPLANE IN 
ORDER TO GAIN/REGAIN AIRSPEED, E.G. TO GAIN SEPARATION 


0316 
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0320 

0321 

0322 

0323 

0324 

0325 

0326 

0327 

0328 

0329 

0330 

0331 

0332 

0333 

0334 

0335 

0336 

0337 

0338 

0339 

0340 

0341 

0342 


NEGGEF=1 

IF ( LOSELB . GT. -5.0/DEGRD) NEGGEE=0 

IF (PHIB .LE. -30 . 0/ DEGFID .OR. PHIB .GE. 30.0/DEGRD) NEGGEE=0 
SI*EEDA=SQRT ( XEDOTA* #2 + YEDOTA* *2+ZED0TA# #2 ) 

IF (VELB .GT. . 90#SPEEDA ) NEGGEE-0 

C 

C 

C — SELECT THE MANEUVER WHICH YIELDS THE SMAL.LEST VARIATION IN ROTB 
C 

I F ( ABS ( ROTBNG-ROTPRV ) .LT. ABS (ROTB-ROTPRV ) 

1 .AND. NEGGEE .EQ. 1) THEN 

GLFVLB= -GLVNEG 
ROTB=ROTBNG 
MANVRB=11 
GOTO 998 
ELSE 

GLFVLB=GLVPOS 
MANVRB=6 
GOTO 998 
END IF 
C 

END IF 

C 21 APR I 

C 

C *** SEE IF WE ARE IN TROUBLE AND NEED AN EVASIVE MANEUVER 
C 
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0343 IF(EPSD.GT. 120. .AND. FLAMBD. GT. 120. ) THEN 

0344 GLEVLB=0. 9 

0345 I F ( CAS . GT . 330 . ) GLEVLB=1 . 

0346 ROTB^ROTB— PIDV2 

0347 MANVRB=8 

0348 END IF 

0349 C 

0350 998 CONTINUE 

0351 WR I TE < # » 491 ) TNOV, MANVRB . ROTB*DEGRD, GLEVLB , RANGE 

0352 C 

'0353 WRITE <77, 49 1 ) TNOW , MANVRB , ROTB*DEGRD» GLEVLB, RANGE 

0354 C 

0355 491 FORMAT < ’ SELCTB S DEC I S I ON* FI 5. 2 . 1 5 , F 10. 2 . F 1 0 . 2 , F 1 2 . 1 . / / ) 

0356 C 

0357 I CNT= I CNT+3 

0358 C 

0359 C »** END OF REGULAR DECISION MAKING PART 

0360 C 

0361 CALL THRTLB 

0362 END IF 

0363 C 

0364 C END OF LEAD/LAG IF 1111111111111111111111111111111111111111111111 

0365 C 

0366 999 RETURN 

0367 C 

0368 C MANEUVER CODES 

0369 C 

0370 C MANVRB=1 DIVE RECOVERY ACTIVE 

0371 C MANVRB=2 DIVE RECOVERY INITIATED 

0372 C MANVRB=3 LEAD PURSUIT 

0373 C MANVRB=4 PURSUIT < FLAMBD= TO : EPSD= TO ) 

0374 C MANVRB=5 LAG PURSUIT 

0375 C MANVRB=6 INTERCEPT TRAJECTORY 

0376 C MANVRB=7 POINTING ALGORITHM 

0377 C MANVRB=8 TAIL-CHASE EVASION 

0373 C MANVRB=9 PURSUIT < FLAMBD® TO ;EPSD= TO ) 

0379 C MANVRB=1 1 NEGATIVE G’S 

0380 END 
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SUBROUTINE OVRSHT ( XEA * YEA , ZEA , XEDOTA * YEDOTA , ZEDOT A * DMTRXA ) 

C 

COMMON/RNDMAN/ MLDECS , MAN I N I * I SLCTR 
C 

COMMON/ CNSTNS/ DT , TBEGN * TNOV * P I , P I DV2 * P I DV4, TVOP I , DEGRD, RADDG, G, 

1 VAR (20) * I VAR (20) , TEND 

C 

COMMON/ CONTRL/ MSTOP , IPRINT 
C 

COMMON/ COMNDB/ I CMNVB* GLEVLB, ROTB, MANVRB 
C 

COMMON/ DATA 1 B/ XEB , YEB , ZEB , XEDOTB * YEDOTB * ZEDOTB * XEDDTB , YEDDTB * 

1 ZEDDTB, PSIB- THETAB, PH IB* UB, VB, WB, PB» QB, RB , A1B, A2B, 

2 A3B * A4B * VELB , VHORB 
C 

COMMON/ DATA2B/ ALFAB * BETAB * CBARB (3,3) , CDB , CLALFB * DMTRXB ( 3* 3) * DRAGB , 

1 L I FTB , LODMXB , LODSTB , MACHB , RHOB * SPECEB , SB , THRSTB , 

2 PSUBSB* TPOSB, I N I ZB , UE I TB * CSB , CLB , PS I BRB , THETBB , 

3 AN 1 B , AN2B , AN3B , MASSB 
C 

REAL L I FTB » LODMXB , LODSTB * MACHB * MASSB 

C 
C 

COMMON/ PO I NTP/ XEA I M , YEA I M , ZEA I M * PCOMB , QCOMB , RCOMB , I POINT 
C 

DIMENSION CMPL (3.3) 

C 

COMMON/ RELVAR/ LOSELA , LOSELB , LOSAZA , LOS A ZB , LOSANA * LOSANB , 

1 LSDOTA * LSDOTB , DEV AN A, DEVANB , DVDOTA , DVDOTB , 

2 RANGE , RRATE , XAINB* YAINB,ZAINB, XBINAfYBI NA , 

3 ZB I NA f ANGOFA f ANGOFB 

C 

REAL LOSELAf LOSELBf LOSAZAf LOSAZBf LOSANA, LOSANBf 

1 LSDOTA, LSDOTB 

C 

C THIS SUBROUTINE GENERATES COMMANDS TO THE B AIRCRAFT TO 

C FORCE AN OVERSHOOT OF THE A AIRCRAFT . THIS MANEUVER COMBINES ROLL 

C IN THE DIRECTION OF A, COMBINED WITH A DIVE 

C 

C THIS MANEUVER SHOULD BE INITIATED UNDER THE FOLLOWING CONDITIONS: 

C RANGE <=3500 FT: EPSD 

C 

C INITIALIZATION SECTION 

C 

I F ( MAN I N I . EQ. -1 ) THEN 
T I MREQ= 15. 

MLDECS= I NT ( T I MREQ/ DT ) 

T I N I T= TNOW 
TQU I T= TNOW+T I MREQ 
MAN I N I =0 
END IF 
C 

DTPRED=0 . 0 

XEXA= XEA + DTPRED#XEDOTA 
YEXA= YEA + DTPRED# YEDOTA 
ZEXA= ZEA + 10000. 

TAXE=XEXA~XEB 
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TAYE=YEXA-YEB 
TA7.E=ZEXA-ZEB 
VH0R2=XED0TB # #2+YED0TB * * 2 
VHORB=SaRT ( VH0R2 ) 

VEI 2=VH0R2+ZED0TB * #2 
VEI .B=SQRT C VEL2 ) 

DZ^(XED0TB*ZED0TB*TAXE+YED0TB#ZED0TB*TAYE~VH0RB**2*TAZE) /VELB 
DY r -YEDOTB*TAXE+XEDOTB*TAYE 
I F ( DZ. EQ. 0. • AMO. DY.EQ.0.) THEN 
ROTB-0. 

ELSE IF COY. EQ. 0. ) THEN 
I F < DZ. GT. 0. ) ROTB=0. 

IFCDZ.LT.0.) ROTB=PI 
ELSE 

R0TB=ATAN2(DY,DZ) 

END IF 
C 

C SELECT MAXIMUM G TRUN 
C 

GLEVLB=0. 95 
TPOSB r 2.0 
ICMNWB=1 
C 

C CHECK FOR MANEUVER TERMINATION CONDITIONS 
C 

I F ( TNOW .GE. TQU IT) THEN 
ISLCTR=0 
END IF 

I F ( RRATE . GT* 0. ) THEN 
ISLCTR=0 
END IF 

RETURN 

END 
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000 1 SUBROUT I NE VTOVSH ( XEA . YEA , ZEA * XEDOTA . YEDOTA . ZEDOTA . DMTRXA ) 

0002 C 

0003 COMMON/ RNDM AN/ MLDECS. MAN 1 N I , I SLCTR 

0004 C 

0005 COMMON/ CNSTNS/ DT * TBEGN . TNOW * P UP! DV2 . P I DV4, TVOP I . DEGRD. RADDG, G. 

0006 1 VAR ( 20 ) . I VAR ( 20 ) * TEND 

0007 C 

0008 COMMON/ CONTRL/ MSTOP t I PR I NT 

0009 C 

0010 COMMON/ COMNDB/ 1CMNWB. GLEVLB. ROTB. MANVRB 

0011 C 

00 1 2 COMMON/ DATA 1 B/ XEB . YEB . ZEB . XEDOTB . YEDOTB * ZEDOTB . XEDDTB , YEDDTB . 

0013 1 ZEDDTB . PS I B . THETAB * PH I B . U8 » VB . WB . PB . QB . RB. A 1 B » A2B . 

0014 2 A3B . A4B > VELB * VHORB 

0015 C 

0016 COMMON/ DATA2B/ ALFAB * BETAB * CBARB (3.3) , CDB . CLALFB. DMTRXB (3*3) . DRAGB * 

0017 1 LI FTB . LODMXB . LODSTB . MACHB . RHOB . SPECEB , SB . THRSTB , 

0018 2 PSUBSB. TPOSB. I N I ZB , WE I TB , CSB , CLB . PS I BRB . THETBB , 

•0019 3 AN1B.AN2B. AN3B * MASSB 

0020 C 

002 1 REAL L I FTB . LODMXB , LODSTB , MACHB . MASSB 

0022 C 

0023 C 

002 4 COMMON/ PO I NTP/ XEA I M * YEA I M , ZEA I M * PCOMB * QCOMB . RCOMB » I PO I NT 

0025 C 

0026 DIMENSION CMPLC3.3) 

0027 C 

0028 COMMON/ RELVAR/ LOSELA. LOSELB , LOSAZA * LOSAZB* LOSANA. LOSANB . 

0029 1 LSDOTA . LSDOTB , DEVANA . DEVANB , DVDOTA , DVDOTB - 

0030 2 RANGE* RRATE. XAINB.YA1NB.ZAINB.XBINA.YBINA. 

0031 3 ZB I NA. ANGOFAf ANGOFB 

0032 C 

3033 REAL LOSELA. LOSELB, LOSAZA* LOSAZB * LOSANA * LOSANB. 

0034 1 LSDOTA. LSDOTB 

0035 C 

0036 C THIS SUBROUTINE GENERATES COMMANDS TO THE B AIRCRAFT TO 

0037 C CLIMB IN THE VERTICAL PLANE FOR A SPECIFIED TIME. 

0038 C THIS CAN BE USED TO FORCE AN OVERSHOOT OF THE A AIRCRAFT 

0039 C WHEN IT IS DIVING ONTO B 

0040 C 

0041 C 

0042 C THIS MANEUVER SHOULD BE INITIATED UNDER THE FOLLOWING CONDITIONS: 

<0043 C RANGE <=3500 FT: EPSD 

0044 C 

0045 C INITIALIZATION SECTION 

0046 C 

0047 I F ( MAN I N I . EQ. -1 ) THEN 

0048 T I MREQ= 15. 

>0049 MLOECS= I NT ( T I MREQ/ DT ) 

0050 T I N I T= TNOW 

0051 TQU I T= TNOV+T I MREQ 

0052 MANINI=0 

0053 END IF 

0054 C 

0055 ROTB=0. 

0056 C 

0057 C SELECT LOW G IN ORDER NOT TO LOSE TOO MUCH ENERGY 
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C 

GLEVLB= 1 - 5/ LODMXB 
TPOSB * 2.0 
I CMNWB= 1 
C 

C CHECK FOR MANEUVER TERMINATION CONDITIONS 
C 

I F < TNOW . GE. TQU I T) THEN 
I SLCTR=0 
END IF 

I F ( RRATE . GT. 0. ) THEN 
I SLCTR=0 
END IF 

RETURN 

END 
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SUBROUT I NE OPSTRN ( XEA * YEA* ZFA * XEDOTA* YEDOTA* ZEDOTA* DMTRXA) 

C 

COMMON/ RNDMAN/ MLDECS * MAN I N I * ISLCTR 
C 

COMMON/ CNSTNS/ DT * TBEGN * TNOW* P I * P I DV2 * P I DV4* TWOP I * DEGRD* RADDG, G» 

1 VAR ( 20 ) * I VAR ( 20 ) * TEND 

C 

COMMON/ CONTHL/ MSTOP * I PR I NT 
C 

COMMON/ COMNUB/ I CMNWB . GLEVLB * ROTB , MANVRB 
C 

COMMON/ DATA 1 B/ XEB, YEB * ZEB * XEDOTB * YEDOTB * ZEDOTB , XEDDTB, YEDDTB* 

1 ZEDDTB* PS I B * THETAB * PH I B * UR * VB * WB* PB * QB * RB * A IB . A2B * 

2 A3B i A4B. VELB * VHORB 

C 

COMMON/ DATA2B/ ALFAB * BETAB * CBARB ( 3 * 3 ) * CDB * CLALFB * DMTRXB ( 3* 3 ) * DRAGB * 

1 L I FTB * LODMXB * LODSTB * MACHB * RHOB * SPECEB * SB f THRSTB * 

2 PSUBSB* TPOSB* I N I ZB * WE I TB * CSB * CLB * PS I BRB * THETBB * 

3 AN 1 B * AN2B , AN3B * MASSB 

C 

REAL L I FTB , LODMXB * LODSTB * MACHB * MASSB 

C 
C 

COMMON/ PO I NTP/ XEA I M * YEA I M * ZFA I M * PCOMB * QCOMB * RCOMB * I PO I NT 
C 

DIMENSION CMPL (3*3) 

C 

COMMOM/RELVAR/LOSELA* LOSELB* LOSAZA* LOSAZB • LOSANA * LOSANB * 

1 LSDOTA ♦ LSDOTB* DEVANA* DEVANB* DVDOTA* DVDOTB, 

2 RANGE * RRATE * XA I NB * YA I NB * ZA I NB , XB I NA * YB I NA * 

3 ZB I NA » ANGOFA* ANGOFB 

C 

REAL LOSELA* LOSELB * LOSAZA* LOSAZB - LOSANA* LOSANB * 

1 LSDOTA* LSDOTB 

C 

C THIS SUBROUTINE GENERATES COMMANDS TO THE B AIRCRAFT TO 

C TURN OPPOSITE TO A FOR A DURATION OF 15 SECONDS 

C OR UNTIL A IS WITHIN A 60 DEG CONE ANGLE 

C IT IS USED TO GENERATE A NOSE-TO-NOSE CONVERSION (B CONVERTS 

C TO A’S NOSE) 

C 

C THIS MANEUVER SHOULD BE INITIATED UNDER THE FOLLOWING CONDITIONS: 

C 

c — . 

C INITIALIZATION SECTION 

C 

I F ( MAN INI. EQ. -1 ) THEN 
T I MREQ= 20. 

MLDECS= I NT ( T I MREQ/ DT ) 

T I N I T= TNOW 
TQUIT= TNOW+T I MREQ 
MAN I N I -0 
C 

c 

DTPRED=0. 0 

XEX A= XEA + DTPRED#XEDOTA 
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YEXA= YEA + DTPRED# YEDOTA 
ZEXA= ZEA 

TAXE=XEXA-XEB 

C 

C NOTE 

C 

TAYE=- ( YEXA-YEB ) 

C 

TAZE=ZEXA-ZEB 
VH0R2=XED0TB*#2+YED0TB * #2 
VHORB=SQRT ( VH0R2 ) 

VE1 .2“VH0R2+ZED0TB**2 
VELB=SQRT ( VEL2 ) 

DZ r <XEDOTB«ZEDOTB*TAXE+YEDOTB*ZEDOTB*TAYE-VHORB#*2#TAZE>/ VELB 
DY - -YEDOTB*TAXE+XEDOTB*TAYE 
I F ( DZ. EQ- 0. -AND. DY.EQ.0.) THEN 
ROTB=0. 

els;e if (dy.eq.0.) then 

IFCDZ.GT.0. > ROTB=0 . 

IF(DZ.LT.0.) ROTB=P I 
ELSIE 

R0TB=ATAN2<DY,DZ> 

END! F 
C 

END IF 
C 

C SELECT MAXIMUM G TURN 
C 

GLEVLP=0 . 95 
TPOSB r 2.0 
CAS=VLLB#0. 5925 
CAS=TAS*SQRT ( RHOB/0 . 0023768 ) 

VCORNH=330. 

IF (CAS . GT. VCORNR) GELVLB-LODSTB/ LODMXB 
I CMNWR= 1 
C 

C CHECK FOR MANEUVER TERMINATION CONDITIONS 
C 

I F ( TNOW . GE. TQU I T ) THEN 
I SLCTR=0 
END IF 
C 

C FLAMBD=LOSANB #DEGRD 

C I F ( RRATE • GT. 0. ) THEN 

C I SLCTR=0 

C END I F 

RETURN 

END 
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SUBROUTINE CL I MB ( XEA • YEA , ZEA, XEDOTA, YEDOTA, ZEDOTA , DMTRXA ) 

C 

COMMON/ RNDMAN/ ML DECS , MAN I N 1 * ISLCTR 
C 

COMMON/ CNSTNS/DTf TBEGN * TNOW, P I , P I DV2* P I DV4 * TWOP I , DEGRD, RADDG » G* 

1 VAR (20) , I VAR (20) * TEND 

C 

COMMON/ CONTRL/ MSTOP * I PR I NT 
C 

COMMON/ COMNDB/ I CMNWB , GLEVLB , ROTB , MANVRB 
C 

COMMON/ DATA 1 B/ XEB , YEB , ZEB * XEDOTB * YE DOT B * ZEDOTB * XEDDTB, YEDDTB* 

1 ZEDDTB * PS I B » THETAB * PH I B » UIU VB » WB* PB f QB * RB , A IB , A2B * 

2 A3B , A4B • VELB * VHORB 

C 

COMMON/ DAT A2B/ALFAB, BETAB * CB ARB (3*3) * CDB * CLALFB , DMTRXB ( 3 , 3) ♦ DRAGB * 

1 L I FTB , LODMXB , LODSTB , MACHB , RHOB , SPECEB , SB , THRSTB , 

2 PSUBSB* TPOSB, IN I ZB, WE I TB * CSB , CLB , PS I BRB , THETBB * 

3 AN 1 B f AN2B * AN3B , MASSB 
C 

RFAL LI FTB, LODMXB, LODSTB , MACHB , MASSB 

C 

COMMON/ PO I NTP/XEA I M, YEA I M , ZEA I M , PCOMB , QCOMB , RCOMB , I PO I NT 

C 

DIMENSION CMPL (3,3) 

C 

COMMON/ RE L VAR/ LOSEL A , LOSELB , LOSAZA * LOSA7.B , LOS AN A , LOSANB * 

1 LSDOTA, LSDOTB, DEVANA, DEVANB, DVDOTA, DVDOTB , 

2 RANGE , RRATE , X A I NB , Y A I NB ,ZAINB, XBINA»YBINA, 

3 ZB I NA, ANGOFA, ANGOFB 
C 

REAL LOSELA, LOSELB, LOSAZA, L0SA7.B , LOSANA , LOSANB, 

1 LSDOTA, LSDOTB 

C 

C THIS SUBROUTINE GENERATES COMMANDS TO THE B AIRCRAFT TO 

C CLIMB IN THF VERTICAL PLANE FOR A SPECIFIED TIME. 

C THIS CAN BE USED TO FORCE AN OVERSHOOT OF THE A AIRCRAFT 

C WHEN IT IS DIVING ONTO B 

C 

c 

C THIS MANEUVER SHOULD BE INITIATED UNDER THE FOLLOWING CONDITIONS: 

C RANGE <=3500 FT; EPSD 

C 

C INITIALIZATION SECTION 

C 

I F ( MAN INI. EG. -1 ) THEN 
T I MREQ=5 . 

MLI)ECS= I NT ( T I MREQ/DT ) 

T I N I T= TNOW 
TQU I T= TNOW+T I MREQ 
MAN I N I =0 
END IF 
C 

ROTB=0. 

C 

C SELECT LOW G IN ORDER NOT TO LOSE TOO MUCH ENERGY 
C ELSE SELECT HIGH G IN ORDER TO GAIN ALTITUDE RAPIDLY 
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C 

GLEVLB=0.95 
TPOSB ? 2.0 
ICMNWB^l 
C 

C CHECK FOR MANEUVER TERMINATION CONDITIONS 
C 

IFCTNOW -GE. TQU I T> THEN 
ISLCTR=0 
END IF 

C I F ( RRATE . GT . 0 . ) THEN 

C I SI CTR=0 

C END IF 

RETURN 

END 


COMMAND QUALIFIERS 

FORTRAN/ L I ST/ SHOW^ NOMAP APPEND I X 

/ CHECK- ( NOBOUNDS, OVERFLOW , NOUNDERFLOW ) 

/ DEQUG= ( NOSYMBOLS, TRACES ACK ) 

/ ST ANDARD- ( NOSYNTAX , NOSOURCE_FORM ) 

/ SHOW= ( NOPRFPROCESSOR* NO I NCLUDE* NOMAP * NOD I CT I ONARY » S I NGLE ) 

/ WARN I NGS= ( GENERAI . , NODECLARAT I ONS , NOULTR I X ) 

/CONTI NUAT I ONS= 1 9 / NOCROSS__REFERENCE / NOD_L I NES / NOEXTEND_SOURCE 

/NOG FLOATING / 1 4 /NOMACHINE CODE /OPTIMIZE 


COMPILATION STATISTICS 

12.01 seconds 
13.04 seconds 
1071 

552 pages 


Run Times 
Elapsed Times 
Page Faults: 
Dynamic Memory: 


/ FY7 


A-15 
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